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PREFACE 


' 


This  report  describes  the  computation  of  range  and  bearing  to 
target  aircraft  from  BCAS  measurements.  Four  different  cases  are 
treated  which  differ  in  respect  to  the  assumptions  made  about  the 
ATCRES  ground  radar  environment.  In  each  case,  only  the  minimal 
configuration,  involving  the  smallest  possible  number  of  radars 
and  targets  necessary  to  obtain  a solution,  is  considered.  All 
computations  are  performed  on  the  basis  of  measurements  at  a single 
point  and  assume  no  a priori  knowledge  of  the  configuration.  All 
the  algorithms  described  and  used  in  the  simulations  were  developed 
at  TSC  and  implemented  in  FORTRAN  on  the  DEC  PDP-10  time-shared 
computer  system.  The  only  exception  is  the  search  algorithm  for 
the  solution  of  the  single-site  configuration,  whose  key  steps  were 
implemented  on  a programmable  hand  calculator.  The  simulations  were 
conducted  using  simulated  random  measurement  variations  consistent 
with  the  measurement  errors  observed  in  flight  tests  of  the 
Litchford  BCAS  feasibility  model. 

The  work  reported  was  sponsored  by  Richard  F.  Bock  of  ARD-250. 
The  draft  was  reviewed  by  Robert  W.  Wisleder,  DTS-531,  and  by  Owen 
E.  Mclntire  and  John  Brennan,  both  ARD-251.  Much  of  the  computer 
programming  was  performed  by  Henrick  0.  Lind,  Barry  I.  Fink,  and 
Ronit  Procaccia  of  Kentron  International,  Ltd.  The  typing  of  the 


draft,  although  complex,  was  well  done  by  Barbara  Davis. 
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1.  INTRODUCTION 


1 . 1 OVERVIEW 

This  report  presents  analyses  and  a design  of  a set  of 
algorithms  for  determining  the  positions  of  intruder  aircraft  and 
ground  interrogators  relative  to  the  BCAS  system  from  the  measure- 
ments that  the  BCAS  system  makes.  Only  the  static  cases  are  con- 
sidered, i.e.,  the  determination  of  positions  on  the  basis  of 
measurements  made  at  the  current  position. 

Five  distinct  cases  are  considered: 

1)  Active  interrogations  of  the  intruder  only. 

2)  Location  of  an  intruder  when  there  is  only  one  ground 
radar,  but  active  interrogation  is  permitted. 

3)  Passive  location  of  the  intruder,  given  two  radars  with 
azimuth  references. 

4)  Passive  location  of  the  intruder,  given  two  radars,  of 
which  only  one  has  azimuth  reference  signals. 

5)  Passive  location  of  the  intruder,  given  two  radars,  of 
which  neither  has  an  azimuth  reference  signal.  (In  both  this  and 
the  previous  case  (4),  two  intruder  aircraft  must  be  repdying  to 
both  radars  or  no  solution  is  possible.) 

For  each  of  these  situations,  only  the  simplest  case  is 
treated,  that  is,  the  case  involving  the  smallest  number  of  air- 
craft necessary  to  solve  for  the  unknown  quantities  describing  the 
configuration.  For  instance,  in  case  (2),  where  the  position  of 
the  intruder  is  determined  by  combining  passive  mode  measurements 
from  one  radar  with  active  mode  measurements,  the  solution 
consists  of  the  range  to  radar  and  bearing  to  only  a single  intruder. 


j 


If  several  targets  are  present,  improved  results  can  be  expected 
because  the  system  of  equations  to  be  solved  then  becomes  over- 
determined. (Clearly,  each  single-intruder  solution  gives  a 
separate  value  for  the  distance  to  the  radar. ) The  computed 
results  can  then  be  expected  to  be  more  accurate  because,  in 
effect,  the  measurement  errors  are  reduced  by  averaging. 

Similar  improvements  in  the  effective  accuracy  of  computing 
intruder  positions  are  to  be  expected  when  aircraft  positions  are 
determined  at  a succession  of  points  and  the  tracks  are  smoothed 
since,  again,  the  errors  are  reduced  by  averaging  over  time.  The 
consideration  of  these  more  complex  cases  is  outside  the  scope  of 
the  present  report. 

1.2  THE  FUNDAMENTAL  EQUATIONS 

The  basic  measurements  performed  by  passive  BCAS  are  of 
(1)  the  differential  time  of  arrival  (TOA),  i.e.,  the  time  dif- 
ference between  the  reception  of  an  SSR  interrogation  and  the 
reception  of  its  elicited  target  reply,  and  of  (2)  the  differ- 
ential azimuth  (DAZ),  i.e.,  the  angle  at  the  interrogator  between 
the  protected  aircraft  and  the  target.  Also,  for  SSR  sites  that 
have  been  suitably  modified,  own  azimuth  (OAZ),  i.e.,  the  magnetic 
beaming  of  own  from  the  radar  site,  can  be  determined.  The 
quantities  measured  are  indicated  in  Figure  1.2-1. 

The  purely  active  BCAS  measures  the  TOA  of  replies  to  its 
own  interrogations.  This  is  the  round-trip  travel  time  of  the 
signals  between  the  BCAS  aircraft  and  the  intruder;  in  effect  it 
amounts  to  a direct  measurement  of  slant  range. 


NORTH 


NORTH 


FIGURE  1.2-1.  SINGLE-SITE  GEOMETRY* 


Both  the  active  and  passive  BCAS  obtain  own  altitude  from 
the  on-board  encoding  altimeter  and  the  intruder  altitude  by 
decoding  the  mode  C replies  received. 

Geometrically,  the  measurements  from  one  SSR  relate  the 
position  of  the  target  aircraft  to  that  of  the  BCAS  in  the  fol- 
lowing way: 

For  a given  separation  d between  the  radar  and  the  BCAS, 
the  TOA  determines  an  ellipsoid  of  revolution  on  which  the  tar- 
get must  lie.  The  radar  and  the  BCAS  are  at  the  foci  of  this 
ellipsoid.  The  differential  azimuth  determines  a vertical  plane 
passing  through  the  radar  and  the  target.  The  target  altitude 
determines  its  horizontal  plane. 
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The  following  equations  describe  the  relationships  among 
the  various  parameters  analytically.  The  subscript  i designates 
one  SSR  out  of  the  available  set. 

The  measured  quantities  for  passive  BCAS  are: 


= OAZ,  azimuth  of  the  BCAS  from  SSR 

= differential  azimuth  (DAZ) 

= the  time  of  arrival  (TOA) 

= altitude  of  target 

= altitude  of  BCAS 

U 

The  (initially  unknown)  quantities  describing  the  configu- 
ration are: 


H 

H_ 


slant  distance  between  SSR  and  target 


= slant  distance  between  SSR  and  BCAS 


R = slant  distance  between  BCAS  and  target 
0 = bearing  of  target  from  BCAS 

Then  the  TOA,  by  definition,  satisfies 


li  - ? <si  + R - Di>  ' 


1.2.1 


where  c is  the  velocity  of  light.  For  convenience,  all  expres- 


sions hereafter  will  be  written  in  terms  of  T^  = c t^  = S^  + R 


- V 


All  angles  are  measured  in  the  horizontal  plane.  By  the 
law  of  cosines,  DAZ  satisfies 


cosct^  = 


? ? ? 

Sf  + DT  - R - 2 HH 
1 1 o 


1.2.2 


2 |/(s2  - H2)(d2  - H2) 


Finally,  it  follows  from  the  law  of  sines  that 


When  active  measurements  are  made,  the  TOA  of  the  replies 
elicited  by  the  active  interrogations  is  tQ  and  satisfies 

t = — . 1.2. U 

o c 

The  sections  that  follow  describe  the  solution  of  these  basic 
equations  for  the  various  cases  described  above.  Each  is  different 
in  that  the  set  of  measurements  - and,  hence,  the  set  of  unknowns 
for  which  one  must  solve  - is  different.  The  solutions  are 


described  qualitatively  in  terms  of  the  locus  curves  determined 
by  the  set  of  measurements  based  on  the  replies  to  a given  radar. 
Then  actual  solution  algorithms  and  simulation  results  are  given. 

1.3  ACTIVE  BCAS 

The  case  of  purely  active  BCAS  will  not  be  treated  in  detail. 
Only  the  range  to  the  intruder,  but  not  the  bearing,  can  be 
obtained.  Flight  test  data  show  the  measurements  to  be  consistently 
good.  An  illustrative  plot  of  range  measured  by  active  interroga- 
tion is  shown  in  Figure  1.3-1. 

The  calculations  of  static  separation  are  trivial.  The  slant 
separation  is 


1.3.1 


1.3.2 


2.  THE  LOCUS  CURVE  APPROACH  TO  DETERMINING  POSITION  OF  OTHER 


The  following  discussion  is  intended  to  describe  qualitatively 
the  nature  of  the  solutions  that  may  be  obtained  to  the  BCAS  equa- 
tions. A sequence  of  essentially  geometric  arguments  will  be  used 
to  construct  the  solutions  by  graphical  means.  In  no  sense  is  it 
proposed  that  graphic  techniques  of  solution  are  appropriate  for  an 
operational  BCAS  system.  However,  the  graphical  solutions  do 
illustrate  the  circumstances  under  which  different  kinds  of  solu- 
tions are  obtained,  indicate  the  error  sensitivity  of  the  solutions 
in  various  geometries,  and  suggest  practical  algorithms  that  might 
be  implemented  to  find  the  solutions. 

2.1  TWO  RADARS  WITH  AZIMUTH  REFERENCE,  THE  PLANAR  CASE 

LOCUS  CURVE  SOLUTION 

2.1.1  Construction  of  the  Locus  Curve 

Assume  that  the  radars  and  both  the  BCAS  aircraft  (OWN)  and 
the  potential  threat  aircraft  (OTHER)  are  constrained  to  a hori- 
zontal plane.  The  measurements  based  on  one  radar  and  the  target 
yield  8,  a,  and  T,  as  indicated  in  Figure  2.1-1. 


FIGURE  2.1-1.  BCAS  MEASUREMENTS  FROM  ONE  RADAR  FOR  ONE 
TARGET.  8,  a,  AND  T ARE  MEASURED,  r,  <J> , d AND  s ARE 
UNKNOWN 
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The  radar-to-OWN  distance  Is  unknown.  If  it  were  known,  the 
position  of  OTHER  could  be  obtained  by  the  following  geometric 
construction: 

With  the  position  of  the  radar  and  the  OWN  aircraft  fixed, 
the  TOA  determines  an  ellipse  on  which  the  target  must  lie.  The 
radar  and  OWN  are  the  foci  of  this  ellipse;  its  major  axis  is  of 
length  T + d.  [This  follows  from  the  familiar  property  of  an 
ellipse:  an  ellipse  is  a curve  such  that  the  sum  of  the  distances 

from  any  point  on  the  curve  to  two  given  points,  called  the  foci, 
is  constant  - i.e.,  such  that  s + r = T + d,  with  T + d given]. 
Furthermore,  the  target  must  lie  in  a direction  from  the  SSR  deter- 
mined by  its  azimuth  measurement,  i.e.,  on  the  line  eminating  from 
the  radar  in  the  direction  (3+a)  from  north.  Thus,  one  may  con- 
struct the  ellipse  and  Intersect  it  with  a line  of  known  direction. 
The  target  must  lie  at  the  point  of  intersection,  assuming  the 
correct  value  of  the  radar-to-own  distance  was  used  in  the  con- 
struction. 

In  fact,  d is  initially  unknown.  Thus  all  that  can  be  deter- 
mined from  the  measurements  of  T,  ct,  and  g is  the  set  of  positions 
for  other  that  is  consistent  with  some  d,  i.e.,  the  locus  of 
possible  target  positions.  A graphical  technique  for  constructing 
the  locus  is  shown  in  Figure  2.1-2.  Two  sets  of  actual  locus 
curves  for  the  planar  case  are  shown  in  Figures  2.4-1  and  2.4-16. 

Several  properties  of  the  locus  curve  are  worth  observing  at 
this  time.  In  generalized  form,  they  will  provide  the  motivation 
for  many  of  the  arguments  In  the  more  complex  cases  to  follow. 
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a,  b,  c,  d are  successive 
assumed  radar  positions 


FIGURE  2.1-2.  GEOMETRIC  CONSTRUCTION  OF  TARGET  LOCUS 
CURVE  FOR  PLANAR  CASE. 

1)  The  locus  curve  is  confined  to  that  sector  of  the  plane 
bounded  by  the  line  of  position  (LOP)  of  OWN  to  radar  and  a line 
intersecting  the  LOP  at  OWN  at  the  differential  azimuth  angle  a. 

Any  point  not  in  this  sector  can  not  be  a target  position  con- 
sistent with  the  given  DAZ. 

2)  Given  any  presumed  target  position  in  the  permitted 
sector,  on  the  locus  curve  or  not,  the  corresponding  radar  position 
can  be  determined  geometrically  by  drawing  a line  through  that 
point  to  intersect  the  LOP  at  the  angle  a.  The  point  of  inter- 
section is  the  radar  position. 

3)  Given  the  assumed  positions  of  the  target  and  the  radar 
relative  to  OWN,  the  corresponding  TOA  can  be  computed.  This 
gives 
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S + R - D . 


This  computed  quantity  Tc  is  equal  to  the  measured  quantity  T (the 
TOA  in  units  of  time  multiplied  by  c,  the  speed  of  light)  If  and 
only  if  the  assumed  target  position  is  on  the  locus  curve. 

k)  The  locus  curve  divides  the  sector  of  the  plane  to  which 
the  target  is  confined  by  the  DAZ  measurement  into  a region  where 
Tc  > T and  a region  where  Tc  < T.  Hence,  by  comparing  T and  Tc  one 
can  determine  on  which  side  of  the  locus  curve  a given  point  is 
without  knowing  the  precise  shape  and  location  of  the  locus  curve 
Itself . 

2.1.2  Finding  the  Target  Position 

When  OWN  is  observing  two  radars  and  the  replies  of  a given 
target  to  both,  then  two  locus  curves  for  the  target  can  be  con- 
structed (as  described  above),  one  corresponding  to  the  set  of 
measurements  from  each  radar.  The  target  must  lie  on  both.  Hence, 
it  must  lie  at  their  intersection. 

Since  the  solution  is  determined  by  the  intersection  of  two 
locus  curves,  it  is  of  interest  to  examine  the  nature  of  the  curves 
and  the  ways  in  which  they  may  intersect. 

Figure  2.1-3a  and  2.1-3b  indicate  schematically  two  kinds  of 
intersections  that  may  occur,  varying  with  respect  to  the  angle  at 
which  the  locus  curves  intersect.  It  is  evident  from  the  diagrams 
that  the  solution  is  more  sensitive  to  measurement  errors  when  the 
curves  intersect  like  those  in  Figure  2.1-3t>  than  in  the  case  for 
the  curves  in  Figure  2.1-3a.  Any  change  in  the  values  of  a,  6, 
and  T used  in  constructing  one  of  the  locus  curves  will  result  in 
a curve  shifted  somewhat  from  the  original.  The  effect  of  such 
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shifts  is  different  in  the  two  cases,  as  shown  in  Figure  2.1-4a 


(a)  (b) 


FIGURE  2.1-3.  TWO  GRAPHICAL  SOLUTIONS  FOR  TARGET  POSITION, 
DETERMINED  FROM  LOCUS  CURVES.  THE  SOLUTION  IN  (a)  IS  LESS 
SENSITIVE  TO  MEASUREMENT  ERROR. 


FIGURE  2.1-4.  SHIFTS  IN  SOLUTION  DUE  TO 
SHIFTS  IN  LOCUS  CURVES. 


The  shift  will  be  relatively  small  when  the  locus  curves  cross 
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approximately  perpendicularly.  The  shift  can  be  very  large  along 

the  direction  of  the  curves  when  the  curves  are  nearly  parallel 

when  they  cross.  ! 

2.2  CONSTRUCTING  LOCUS  CURVES  FOR  THE  GENERAL  THREE-DIMENSIONAL 

CASE 

The  locus  of  possible  target  positions  determined  by  the 
measurement  set  from  a single  ground  radar  can  be  determined  by  a 
similar  geometric  construction  when  the  aircraft  altitudes  are 
nonzero.  However,  the  locus  curves  for  the  general  case  may  have 
qualitatively  different  properties,  which  must  be  understood  if 
some  phenomena  that  may  occur  in  the  general  case  are  to  be 
explained. 

The  steps  in  the  geometric  construction  of  the  locus  curve 
for  a given  combination  of  TOA,  OAZ,  DAZ,  own  altitude  Hq  and 
target  altitude  H are  the  following: 

Again  a succession  of  horizontal  distances  of  OWN-to-radar , 
designated  d,  are  assumed.  The  altitude  of  the  ground  radar  is 
taken  to  be  zero  throughout.  The  geometric  construction  of  the 
target  position(s)  corresponding  to  a particular  OWN-to-radar 
distance  d is  shown  in  Figure  2.2-1.  The  steps  are  the  following: 

1)  With  OWN  Bearing  6,  own  height  Hq,  and  TOA  T known  and 
d assumed,  an  ellipsoid  of  revolution  is  determined  on  which  the 
target  must  lie.  OWN  and  the  assumed  radar  position  are  its  foci. 

The  length  of  its  major  axis  is 

T + D = T + V d2  + H2  . 
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2)  This  ellipsoid  is  assumed  to  intersect  the  horizontal 
plane  determined  by  the  altitude  H of  the  target.  In  the  plane, 
the  intersection  is  an  ellipse  whose  major  axis  lies  along  the 
LOP.  (It  should  be  noted  that  there  will  always  be  some  value  of 
H > Hq  for  which  the  intersection  will  not  exist  - i.e.,  such  that 
the  whole  ellipsoid  is  below  the  target  altitude  plane.) 

3)  In  the  target  altitude  plane,  a line  (the  "DAZ  line")  now 
can  be  drawn  from  the  projected  radar  position  at  an  angle 
determined  by  the  DAZ,  as  before.  It  is  at  this  stage  that  the 
qualitative  differences  between  the  planar  case  and  this  three- 
dimensional  case  become  evident.  In  the  construction  for  the  planar 
case,  the  radar  position  was  inside  the  ellipse  (at  one  focus)  and 
the  line  determined  by  the  DAZ  always  intersected  the  ellipse  at 
exactly  one  point.  In  the  general  three-dimensional  case,  any  of 
the  following  cases  may  occur,  depending  on  the  relationships 
among  the  parameter  values: 

(a)  The  radar  position  projected  vertically  may  fall  inside 
the  ellipse  in  the  altitude  plane  (though  it  will  not  be  at  its 
focus).  Then  the  "DAZ  line"  will  intersect  the  ellipse  at  exactly 
one  point,  as  in  the  planar  case.  At  least  for  values  of  d in  that 
neighborhood,  the  target  locus  curve  will  have  one  branch. 

(b)  The  projected  radar  position  may  fall  outside  the 
ellipse  in  the  altitude  plane  (cf.  Figure  2.2-1).  Then,  depending 
on  the  DAZ,  there  may  be  two  intersections  of  the  ellipse  with  the 
"DAZ  Line",  one  intersection  (at  the  point  of  tangency),  or  no 
intersection  at  all  (when  the  DAZ  is  too  large). 
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The  locus  curves  for  the  general  case  are  the  sequences  of 
points  generated  for  each  value  of  d,  as  d varies  from  zero  to 
infinity.  The  nature  of  the  curves  depends  on  the  relative  values 
of  the  parameters.  Figure  2.2-2  shows  the  cases  that  may  occur. 

The  curve  in  Figure  2.2-2a  is  most  nearly  like  that  obtained 
with  all  altitudes  zero.  The  only  qualitative  difference  is  that 
the  locus  curves  obtained  in  the  planar  case  are  always  concave 
toward  OWN.  Here  the  curve  is  shown  with  an  inflection  away  from 
own  in  the  region  corresponding  to  small  values  of  d (cf.  Figures 
2.4-1  to  2.4-24).  It  will  be  seen  that  this  may  result  in  double 
solutions  for  some  geometries,  which  can  not  occur  in  the  planar  case. 

The  locus  curve  in  Figure  2.2-2b  has  two  branches,  which  come 
about  due  to  the  radar  being  outside  the  ellipse  in  the  altitude 
plane  (cf.  Figure  2.2-1).  The  inner  branch,  i.e.,  the  branch 
nearer  the  LOP,  eventually  joins  the  LOP.  This  occurs  for  a value 
of  d such  that  the  ellipsoid  becomes  so  steeply  inclined  that  the 
projected  radar  position  touches  the  ellipse  in  the  target  alti- 
tude plane  (cf.  Figure  2.2-1).  For  smaller  values  of  d,  the  pro- 
jected radar  position  is  inside  this  ellipse,  and  the  "outer" 
branch  of  the  locus  curve  is  qualitatively  like  the  locus  curve 
in  Figure  2.2-2a. 

It  is  shown  in  the  Appendix  that  a necessary  condition  for  the 
locus  curves  to  have  the  character  of  those  in  Figure  2.2-2b  is 
that  T < H,  the  target  altitude.  Furthermore,  it  is  shown  that 
both  branches  of  the  locus  curve  become  asymptotically  parallel 
to  the  LOP.  Further,  it  is  shown  that,  if  T,  H and  Hq  are  held 
constant  and  the  DAZ  a is  increased,  the  "inner"  branch  of  the  locus 
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curve  moves  away  from  the  LOP  and  the  "outer"  branch  approaches 
the  LOP. 

If  a becomes  sufficiently  large,  the  two  branches  of  the  locus 

/ I 1 

curve  do  not  go  to  infinity,  but  instead  merge,  as  in  Figure  2.2-2c 

(cf.  Figure  2.4-14).  It  is  shown  in  the  Appendix  that  a sufficient 
condition  for  this  to  occur  is  that  T < H sin  a. 

Figure  2.2-2d  shows  a case  when  the  two  branches  of  the  locus 
curve  exist  separately  far  from  OWN,  but  which  then  join,  instead 
of  terminating  separately.  In  terms  of  the  geometrical  process 
lor  constructing  locus  curves  described  previously,  such  a locus 
can  come  about  when  the  following  happens:  Assume  that  the  target 
altitude  is  above  that  of  OWN  by  so  much  that  when  the  radar  is 
assumed  directly  under  OWN,  the  ellipsoid  does  not  reach  the 
target  altitude  plane.  This  will  occur  if  and  only  if 

T < 2 (H-H  ) . 

It  is  shown  in  the  Appendix  that  for  any  T > 0 , as  the  radar  moves 
away  from  OWN,  i.e.,  as  d Increases,  the  ellipsoid  eventually  will 
touch  the  target  altitude  plane  at  a point  over  the  LOP  and  the 
Intersection  will  expand  into  an  ellipse  as  d increases  further. 

For  sufficiently  small  a,  this  ellipse  will  be  cut  by  the  line  from 
the  radar  at  the  angle  determined  by  the  DAZ.  As  d increases,  this 
leads  to  the  same  asymptotic  behavior  as  previously. 

Geometric  intuition  implies  that  a locus  curve  In  the  form  of 
a closed  contour  can  not  occur.  If  a is  so  large  that  no  solution 
exists  for  d -*■  then  there  will  not  be  a combination  of  circum- 
stances such  that  the  ellipse  that  forms  in  the  altitude  plane 
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with  the  radar  position  outside  it  can  first  expand  with  d fast 
enough  to  be  cut  by  the  DAZ  line  and  then  fail  to  expand  fast 
enough,  so  that  ultimately  the  DAZ  line  no  longer  intersects  it. 

2.3  THE  LOCUS  CURVES  AS  BOUNDARIES  OP  REGIONS 

In  the  discussion  so  far,  the  locus  curves  have  been  described 
as  geometrical  constructions  which  can  be  drawn  for  given  sets  of 
measurement  values  T,  a,  £,  H and  H . It  will  be  convenient  in 
what  follows  to  regard  them  from  an  alternative  point  of  view  as 
boundaries  separating  distinct  subregions  of  the  region  where  the 
other  aircraft  may  be  located. 

If  own  azimuth  (that  is,  the  direction  of  the  LOP)  and  the 
differential  azimuth  a are  given,  then  the  potential  locations  of 
the  target  aircraft  are  constrained  to  a region  bounded  by  the  LOP 
and  a line  eminating  from  OWN  at  an  angle  a to  the  LOP.  This  will 
be  referred  to  as  the  solution  region.  It  is  the  region  swept 
out  by  the  line  from  the  radar  at  the  angle  determined  by  (OAZ  + 

DAZ)  as  the  (horizontal)  distance  of  the  radar  to  OWN  ranges  from 
zero  to  infinity. 

The  target  may  be  assumed  to  be  at  any  point  within  this 
region  defined  by  its  polar  coordinates  C r , ) where  is  in  the 
appropriate  interval.  Then  the  radar  location  consistent  with  the 
assumed  target  position  and  the  given  DAZ  can  be  constructed  geome- 
trically by  drawing  a line  through  the  point  (r,4>)  at  an  angle  a to 
the  LOP.  The  radar  position  will  be  at  the  intersection  of  this 
line  and  the  LOP. 

Given  the  position  of  the  target  and  the  radar  relative  to  OWN, 
the  TOA  corresponding  to  the  configuration  can  be  computed  as 
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, (r,4>)  = 


V 


sin/  r2  + H2  + 
sin  a 


r2  + (H-Hq)2 


i 


sin  ( <*+a ) 
sin  a 


r2  + H2 
o 


This  is  seen  to  be  a function  of  r and  <t>  that  is  defined  and  con- 
tinuous everywhere,  and  in  particular  in  the  solution  region, 
within  which  it  has  physical  and  intuitive  meaning. 

At  every  point  on  a target  locus  curve,  Tc  is  equal  to  T,  the 

observed  TCA:  T > T on  one  side  of  the  curve  and  T < T on  the 
c c 

other.  Thus  the  locus  curves  separate  the  solution  region  into 

regions  where  T > T and  T < T.  This  implies  that  the  curves 
c c 

must  either  end  on  a boundary  of  the  solution  region  or  at  infin- 
ity - otherwise  there  would  be  a path  within  the  solution  region 
such  that  along  the  path  the  continuous  finite  function  Tq  - T 
changed  sign  without  ever  going  through  zero. 

In  terms  of  the  process  for  constructing  locus  curves 
described  previously,  it  is  seen  that  Tc  < T on  the  side  of  the 
locus  line  lying  on  the  inside  of  the  constant  - TOA  ellipse  and 
7„  > T on  the  side  corresponding  to  the  outside  of  the  ellipse. 

The  greatest  significance  of  the  observation  that  the  locus 
curves  are  boundaries  lies  is  its  potential  use  for  finding  target 
aircraft  positions  by  search  in  the  solution  region.  An  algorithm 
applicable  to  the  single-site  geometry  is  described  fully  in 
Section  3. A.  An  algorithm  similar  in  principle,  though  more 
complex,  could  be  constructed  for  the  case  of  passive  operation 
with  two  radars  with  azimuth  reference  signals.  The  cases  to  be 
considered  are  identified  and  described  in  Section  H. 
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2.4  THE  NATURE  OF  THE  LOCUS  CURVES 

Plots  of  families  of  locus  curves  are  shown  In  Figures  2.4-1 
to  2.  4-24.  The  shape  of  each  curve  Is  determined  by  the  three 

variables  a,  H/T,  Hq/T.  The  scale  is  determined  by  T.  Only  a 
sampling  of  the  possible  combinations  of  parameters  have  been  used 
In  plotting  the  curves.  The  families  of  curves  to  be  plotted  on 
any  one  graph  were  selected  to  show  the  sensitivity  of  the  curves 
to  changes  in  the  individual  parameters.  The  following  conclu- 
sions appear  to  be  generally  valid: 

1)  For  T > H and  T > Hq,  the  curves  change  relatively  little 
with  changes  in  the  differential  azimuth  (Figures  2-4.1  to  2.4-5). 
Hence,  they  are  insensitive  to  DAZ  measurement  errors. 

2)  Over  a very  wide  range  of  altitudes,  the  curves  change 
very  little  (except  in  the  quadrant  opposite  the  radar)  if  both 
own  and  target  altitudes  change  together,  or,  what  is  equivalent, 
the  radar  altitude  deviates  from  its  nominal  (zero)  value  (Figures 
2.4-6  to  2.4-12).  The  sensitivity  is  least  when  the  target  and  OWN 
are  at  coaltitude.  Hence,  knowledge  of  true  radar  altitude  is  not 
important  to  BCAS  for  determining  target  positions. 

3)  The  curves  do  not  vary  drastically  when  the  relative  alti- 


tudes of  OWN  and  other  change.  The  curves  tend  to  come  closer  to 
the  origin  - i.e.,  to  OWN  when  the  altitude  separation  increases. 
It  is  important  to  note  that  this  happens  whether  OWN  is  above  or 
below  the  target  (Figures  2.4-13  to  2.4-22).  Hence,  it  does  not 
appear  feasible  to  deduce  the  altitudes  relative  to  own  of  other 
aircraft  not  equipped  with  mode  C transponders  from  other  BCAS 


measurements.  It  does  appear  feasible  to  estimate  the  range  and 
bearing  of  such  aircraft  by  assuming  them  to  be  at  coaltitude,  but 
the  horizontal  range  estimates  will  be  somewhat  too  high. 

lj)  The  locus  curves  do  change  quite  drastically  as  the 
various  parameters  change  when  the  TOA  is  small  relative  to  the 
altitude  (Figures  2.4-23  to  2.4-24).  Hence,  solutions  based  on 
measurement  sets  including  small  TOA's  are  likely  to  be  highly 
error-sensitive . 
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FIGURE  2.4-2.  HEIGHT  OF  OWN  = 5000,  HEIGHT  OF  OTHER  = 7000, T = 5-000  MILES. 
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FIGURE  2.4-3.  HEIGHT  OF  OWN  = 10000,  HEIGHT  OF  OTHER  = 11000, T = 3.000  MILES 


FIGURE  2.4-4.  HEIGHT  OF  OWN  = 10000,  HEIGHT  OF  OTHER  =12000, T = 10.000  MILES 


FIGURE  2.11-5.  HEIGHT  OF  OWN  - 28000,  HEIGHT  OF  OTHER  = 31000, T = 5.000  MILES 
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FIGURE  2.4-6.  HEIGHT  OF  OTHER  = 12000,  T = 5.000  MILES 


FIGURE  2.4-7.  HEIGHT  OF  OTHER  = 12000,  T = 5.000  MILES, DAZ  = 5.000 
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FIGURE  2.4-11.  Hi  - HO  = 0,  T = 5.000  MILES, DAZ  = 5-000 


FIGURE  2.4-12.  HI  - HO  = -5000,  T = 5.000 


FIGURE  2.4-13.  HEIGHT  OF  OWN  = 20000,  HEIGHT  OF  OTHER  = 17000, T = 2.000  MILES. 


FIGURE  2.4-111.  HEIGHT  OF  OWN  = 14000,  HEIGHT  OF  OTHER  = 13000, T = 0.200  MILES 


FIGURE  2.4-15.  HEIGHT  OF  OWN  = 25000,  HEIGHT  OF  OTHER  = 26000, T = 0.200  MILES 


FIGIRE  2.4-16.  HEIGHT  OF  OWN  = 0.  HEIGHT  OF 
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FIGURE  2.14-21.  HEIGHT  OF  OWN  = 10000,  HEIGHT  OF  OTHER  = 11000, DAZ 
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3.  SINGLE-SITE  SOLUTION 


3.1  OVERVIEW 

The  single-site  solution  rests  on  measurements  made  upon  the 
replies  to  one  /round  SSR  whose  direction  from  OWN  is  known  and 
upon  the  target  replies  to  active  interrogations  by  OWN.  The 
measurements  available  from  passive  "listening  in"  on  the  ground 
SSR  are 

Own  azimuth  6 
Differential  azimuth  a 
Passive  TCA , where  T = TOA'C  miles 
OWN  altitude  H 

o 

OTHER'S  altitude  H 

From,  this  set  of  measurements,  the  corresponding  single  radar 
locus  curve  can  be  determined.  The  nature  of  this  curve  can  be 
of  any  one  of  the  varieties  described. 

Active  interrogation  yields  a measurement  of  slant  range  to 
OTHER,  R.  With  H and  Hq  known,  the  horizontal  range  to  OTHER  can 
be  determined  as 

r = j R2  - (H-Hq)2  . 

Thus  the  locus  of  possible  positions  of  OTHER  determined  by  active 
interrogations  is  a circle  of  radius  r centered  on  the  vertical 
projection  of  OWN.  The  other  aircraft  must  therefore  be  at  a 
point  of  intersection  of  this  circle  with  the  passive  radar  locus 

curve . 

The  different  possible  kinds  of  interactions  are  shown  in 
Figures  3.1-1  through  3.1-5. 
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Locus  determined  by 

passive  measurements 

of  TOA  = c T 

T > T and  DAZ  = a 


T < T 
c 


Solution 


Locus  determined 
by  active 
' measurement 
of  range  to 
other  = r. 


FIGURE  3.1-1.  NATURE  OF  SOLUTION  FOR  SINGLE-SITE  GEOMETRY 
WHEN  ALTITUDES  ARE  ZERO  (a  < 0). 


T > T 
c 


SSR  _ Tc  > 


Solution  1 


FIGURE  3.1-2.  NATURE  OF  SOLUTIONS  FOR  SINGLE-SITE  GEOMETRY 
WHEN  LOCUS  CURVE  DETERMINED  BY  PASSIVE  MEASUREMENTS  HAS  TWO 
BRANCHES  (a  < 0). 
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Solution  1 


OWN 


FIGURE  3.1-3.  SPECIAL  CASE  OF  SOLUTION,  RELATED  TO  THAT 
IN  FIGURE  3.1-2  (CF.  FIGURE  2. 4-15 X 


FIGURE  3.1-4.  SPECIAL  CASE  OF  SOLUTION,  RELATED  TO  THAT 
IN  FIGURE  3-1-2  (CF.  FIGURE  2.4-14). 
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FIGURE  3.1-5-  CONDITION  WHEN  TWO  SOLUTIONS  CAN  EXIST  DUE  TO 
CIRCLE  DETERMINED  BY  ACTIVE  INTERROGATION  INTERSECTING  OUTER 
BRANCH  OF  PASSIVE  MEASUREMENT  LOCUS  CURVE  TWICE  (a  < 0). 


Two  methods  of  solution  are  possible  for  this  configuration. 

The  first  is  an  algebraic  method  that  analytically  finds  all  roots 
of  the  equation  relating  measurements  to  the  target  position.  The 
equation  is  derived  in  Section  3.2.  This  algorithm  has  been  imple- 
mented as  a FORTRAN  program  written  for  the  PDP-10  time-sharing 
system  at  TSC  and  has  been  successfully  tested. 

The  second  method  of  solution  is  essentially  a search  algorithm 
based  on  the  properties  of  the  locus  curves.  The  algorithm  is  out- 
lined in  Section  3.1*.  The  search  procedure  has  been  implemented  on 
a programmable  hand  calculator  and  successfully  demonstrated.  Due 
to  the  constraints  of  the  calculator  (100  program  steps,  10  memo- 
ries), only  the  search  procedure,  not  the  logic  for  determining 
which  case  pertains,  has  been  implemented  as  a program. 


3.2  PLANAR  CASE,  ALGEBRAIC  SOLUTION 


The  geometry  in  the  single-site  configuration,  with  all 
altitudes  zero,  is  the  simplest  of  the  cases  to  be  considered. 

The  solution  will  be  carried  out  algebraically  to  illustrate  what 
must  be  considered  in  solving  the  BCAS  equations  algebraically. 
The  geometric  configuration  is  defined  by  Figure  3.2-1. 


T = s + r - d 


FIGURE  3.2-1.  SINGLE-SITE,  PLANAR  GEOMETRY 


The  quantities  T,  a and  r are  known,  <j>  and  d are  to  be  found. 
One  may  not  assume  that  d > r.  The  TOA  equation  gives 


T = s + r - d 


3.2.1 


By  the  law  of  sines 


sin  d>  sin  a 


3.2.2 


d _ d _ r 

sin  (l80-a-<}>)  sin  (a+cjj)  sin  a 


3.2.3 


T = il^  sin  + r - iI^  sin  (0+a) 


3.2.4 


m_  vi 

— — sin  a = sin  <J>  - sin  (<}>+a) 


3.2.5 


Now  define 


Then 


Hence 


sin  a = sin  (y~^  ~ sln  (y+|) 
= -2  cos  y sin  ~ 


cos  y 


(r-T)  sin  a 
2 r sin  | 

r-T  « 
— cos  J 


3.2.6 


3.2.7 


3.2.8 


and 


Y 


+ cos"1  cos 

+ cos'1^—  cos 


3-2.9 

3.2.10 


Thus  it  is  seen  that  the  equation  by  which  <p  is  determined  has  two 
roots.  It  remains  to  be  established  which  of  these  roots  is  an 
actual  solution  of  the  geometric  problem. 

Let 


CL  , 

-i 

/ r-T 

= " 2 + 

cos 

(— 

cos 

5) 

3.2.11a 

01 

-1 

( r-T 

a \ 

2 

cos 

\ — 

cos 

2 / 

3.2.11b 

It  is  apparent  that  <J>2  results  in  a position  for  OTHER  on  the  wrong 
side  of  the  LOP  - i.e.,  inconsistent  with  the  sign  of  a.  Hence 
only  is  a valid  solution  to  the  problem  posed.  Thus  this 
problem  is  seen  to  have  a unique  solution,  assuming  that  a solution 
exists.  It  is  to  be  noted  that  there  obviously  exist  sets  of 
values  of  r,  T,  and  a such  that 
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a 

2 


r-T 

r 


cos 


> 1 


Then  there  is  no  value  of  $ that  satisfies  Equation  3.2.5. 

The  fact  that  there  may  be  sets  of  measured  values  for  which 
no  solution  exists  is  important  to  recognize.  Such  a situation 
could  not  arise  if  all  measurements  were  noise-free.  However, 
given  that  the  measurements  are  corrupted  by  independent  random 
noise,  the  possibility  of  "unrealizable"  measurement  sets  exists 
and  must  be  considered  in  designing  any  algorithm  for  determining 
the  position  of  other,  lest  there  be  a program  "crash"  when  the 
situation  occurs. 


3.3  THREE-DIMENSIONAL  CASE  - ALGEBRAIC  SOLUTION 

The  algebraic  equation  to  be  solved  is  derived  as  follows: 

Given 

R = ■yl  r2  + (H-Hq)2  3-3.1 

the  time-of-arrival  equation  expressed  in  range  dimensions  can  be 
written 


T - R +^d2  + H2  =i^s2  + H2 


3-3.2 


which  squares  to 


(T-R)2  + d2  + H2  + 2 (T-R) y d2  + 


2 ? P 

= s^  + Rd  3.3.3 


so  that 


(T-R)2  + d2  - s2  + H2  - H2  = - 2(T-R)Vd2  + H2  3.3.4 

o wo 


By  the  law  of  sines  (cf.  Figure  3-3-1) 

r sin(9-B-ot) 
sin  a 


3.3.5 


and 


r sln(e-B) 

sin  a 


3.3.6 


52 


FIGURE  3.3-1.  SINGLE-SITE  GEOMETRY/PLAN  VIEW. 
Subtraction  of  3. 3-6  from  3*3*5  &nd  trigonometric  manipula 
tion  yields 


d - s = 


sin 


— £sin(9-8-ot)  - sin(e-B)] 


Similarly 


2— - cos(9-B-§)  sin  y 
sin  a 2d 


d + s = ~?r—  sin(e-B-f)  cos(-“-) 
sin  a d d 


3*3*7 


3*3*8 


Mult iolicat ion  of  3.3.7  and  3-3.8  yields 


? 2 Hr2  sin(9-6-a/2)  cos ( 9- 6-a/2 ) sln(a/2)  cos(a/2) 

d - s = - p 

sin  a 

3.3.9 

Let  Y = 26  - 29  + a 
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Then  Equation  3-3-9  can  be  written  as 


,2  _ 2 = r sin  y 

sin  a 


3.3.10 


Substitution  of  Equation  3- 3- 10  in  Equation  3.3-1*  gives 


(T-R)2  + - — -f— - * + H2  - H2 
sin  a o 

so  that,  after  squaring  both  sides 


2 (T-R)  ^d2  + H2  3.3.11 


(T-R ) ^ pn  y+  (H2-H2)2  + ~2  Sin  Y + 2(H2-H2)  (T-R)2 


2 

sin  a 


Now 


+ 2 (H2-H2 ) r2  ^ -Y  = 4 (T-R) 2 (- — Sy”  (e~6~a)  + h2 


(! 


, 2 
sin  a 


3-3.12 


sin2  ( 0-6-a ) = I -_.?0S(2(8-B-«)) 

1 - cos  (-2(8-B-ot)) 
2 

1 - cos  (y+a) 

2 


3.3.13 


1 - cos  y cos  a + sin  y sin  a 
2 


Substitution  of  3.3.13  in  3.3.12,  rearrangement  and  multiplication 
by  sin2  « gives 

sin2y  r**  + siny  ^2r2(H2-H2)  sina|+  cosy|2(T-R)2  r2  cosa| 

+ | (T-R)^  sin2o  + (H2-H2)2  sin2a  + 2 (H2-H2 ) (T-R ) 2 

- 4 (T-R ) 2 H2  sin2a  - 2(T-R)2  r2  } = 0 3.3.14 

This  is  of  the  form 

2 

a^  sin  y + a^  siny  + a^  cosy  + a^  = 0 3.3.15 

2 2 

Transpose  a^  cosy,  square  and  substitute  1 - sin  y for  cos  y. 


54 


The  result  Is 


b^  sin\+  b^  sln^y  + sin2y  + b^  sin  y + bQ  = 0 3-3.16 


where 


b3  = Ur6  (H2-H2)  sin  a 

h t o 0 0 0 ii  o h 0 

b2  = 2r  li(KQ-K  ) sin  a + 2(T-R)  cos  a + (T-R)  sin  a 

+ 2(T-R)2^H2-H2-2H2  sin2a-r2^ 

b]_  = 4r2(H2-H2)  sin  a|(T-R)^  sin2a  + (H2-H2)2  sin2a 

+ 2(T-R)2^H2-H2-2H2  sin2a-r2)} 

bQ  = |(T-R)^  sin2a  + (H2-H2)2  sin2o 

+ 2(T-R)2 (h2-H2-2H2  sin2ot-r2)}2  - ^r4  (T-R) ^ cos2ct 

Equation  3.3.16  is  a quartic  equation  which  may  be  solved  closed 
form  for  sin  y. 

There  may  be  up  to  four  real  roots  of  3.3.16  (i.e.,  up  to 
four  possible  values  of  sin  y)  between  -1  and  1.  There  is  at  most 
one  value  of  9 corresponding  to  each  possible  value  of  sin  y. 

Each  root  is  treated  as  follows: 

Equation  3*3.15  is  used  to  compute  cos  y.  Since  sin  y and 
cos  y are  now  known,  y may  be  computed  uniquely. 

y = 26  - 20  + a so 

9 = \ (2  6+ot-y)  3.3.17 

The  difficulty  now  arises  that  3.3.17  only  determines  4>  up  to  + ir. 
(An  angle  has  two  half-angles.)  However,  there  is  only  one  solu- 
tion 9 consistent  with  the  requirement  that  9 lie  in  the  interval 
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determined  by  8 and  a (cf.  Section  2.1  and  Figure  2.1-1).  For 
1 0—  8 1 < 180° , sgn  ($-6)  = sgn  a.  Thus  the  bearing  angle  0 corre- 

sponding to  every  real  root  of  Equation  3 • 3 - 15  in  the  interval 
(-1,1)  can  be  determined  uniquely.  Then  Equation  3-3.5  is  used  to 
solve  for  d and  3-3-6  for  s. 

Equations  3.3-2  and  3-3-11  are  squared  en  route  to  producing 
3.3-15.  It  is,  then,  possible  to  introduce  extraneous  roots  of 
3.3.16  by  these  steps.  If  the  signs  of  the  left  and  right  hand 
sides  of  3.3.2  or  3.3.3  disagree,  the  root  is  discarded. 

The  fourth-order  polynomial  equation  3-3-15  has  a single  root 
of  multiplicity  four  when  T = R.  The  root,  obtained  by  factoring. 


F2  - H2 

sin  y = 5 — - Sin  a 


3-3-18 


Numerical  techniques  of  solving  the  equation  will  encounter  dif- 
ficulties (due  to  numerical  roundoff  effects)  when  this  condition 
occurs  or  nearly  occurs,  i.e.,  when  T = R. 

The  algorithm  as  coded  tests  for  this  condition,  and  when  it 
occurs,  uses  the  expression  in  3-3-18  as  an  initialization  to 
solve  by  Newton's  method  the  nonlinear  equation  obtained  by  sub- 
stituting in  3-3-2  the  expressions  for  s and  d in  3-3-6  and  3-3-7- 
Alternatively,  the  search  algorithm  described  in  the  next 
section  could  be  used. 

3.4  SINGLE-SITE  GEOMETRY  - SOLUTION  BY  SEARCH 

The  position  of  the  target  aircraft  and  the  distance  to  the 
ground  SSR's  can  be  determined  by  a one-dimensional  search,  based 


on  the  known  properties  that  the  solution  may  have. 


The  cases  are  distinguished  by  the  behavior  of  the  quantity 
Tc  along  the  boundaries  of  the  solution  region.  In  particular, 

Tc  has  to  be  evaluated  at  two  points  to  distinguish  the  differences 
among  the  cases  that  require  different  search  algorithms.  These 
points  are  at 

r = |/r2  - (H-Hq)2  3.4.1 

and  at  <f>  = 0 and  <t>  = 180  - a;  when  <)>  is  defined  by  Figure  3.3-1. 

Tc(r,0)  = T1  = H + R - ^r2  + H2  3-4.2 

and 

Tc(r,l80-a)  = T2  = yr2  + H2  + R - Hq  3-4.3 

The  cases  of  interest  are  defined  in  Table  3-4-1 


TABLE  3-4-1.  DEFINITION  OF  CASES  FOR  SEARCH  LOGIC 
FOR  SINGLE-SITE  GEOMETRY 


T1 

T2 

Case 

Corresponding  Figures 

< T 

> T 

I 

3-1-1 

> T 

> T 

II 

3-1-2,  3-1-3,  3-1-4 

< T 

< T 

III 

3.1-5 

> T 

< T 

Not  possible 

The  condition  that  T^  > T and  T^  < T can  not  occur,  since  this 
would  require  that 

T1  > T > T2  3-4.4 


or 

H + H0  > Vr2  + H2  + ^/r2.  + H2  3.4.5 

which  can  not  occur  for  real  values  of  r. 

I 

Once  the  pertinent  case  has  been  identified,  the  solution 
points  can  be  found  by  search.  The  search  is  based  on  the  observa- 
tion that  Tc  is  greater  than  T on  one  side  of  the  locus  line  and 
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less  than  T on  the  other.  If  (T  -T)  Is  of  opposite  sign  at  two 
points  In  the  solution  space,  the  locus  line  must  lie  between  them. 
The  search  then  proceeds  very  simply.  It  is  known  that  the  solu- 
tion point (s)  must  lie  on  a circle  of  radius  r in  a sector 

0 < <p  < 180°  - a 3-4.6 

for  a > 0 and  a sector 

0 > <(>  > - 180°  + a 3.4.7 

for  a < 0. 

The  search  in  principle  proceeds  by  evaluating  (T  -T)  at  a 
succession  of  values  of  <f>  at  increments  of  A4>.  When  (T  -T)  changes 
sign,  it  is  known  that  the  solution  point  is  within  the  last  pre- 
vious increment.  Accordingly,  the  direction  of  search  is 
reversed,  using  increments  of  <j>  half  as  large  as  before.  As  A<J> 
is  successively  halved,  the  successive  values  of  the  test  points 
approach  the  true  solution  point  and 


FIGURE  3.4-1.  RELATIONSHIP  OF  VARIABLES  USED  IN  SEARCH 
ALGORITHM  FOR  SINGLE-SITE  GEOMETRY  SOLUTION.  CASE  WITH  a 
AND  <(>  POSITIVE  IS  SHOWN. 


^ ^true  ^test^  < I ^ A ^ I 

Thus,  the  Iteration  can  be  terminated  whenever  A<p  is  small 
enough  to  guarantee  the  desired  accuracy. 

Simplified  flow  charts  for  the  solution  by  search  are  shown 
In  Figures  3.^~2  and  3-^-3-  The  search  algorithm  necessarily 
always  converges  when  there  is  exactly  one  solution  point  in  the 
range  being  searched  and  that  solution  arises  because  the  passive 
radar  locus  curve  intersects  the  circle  determined  by  the  active 
radar . 


When  there  is  one  solution  point  due  to  the  passive  radar 
locus  curve  being  tangent  to  the  active  radar  circle,  then  the 
algorithm  will  find  the  solution  only  in  the  unlikely  event  that 
one  of  the  test  points  selected  happens  to  be  the  exact  solution 
point.  Since  the  difference  (T-Tc)  does  not  change  sign  along  the 
circle  on  either  side  of  the  solution  point,  there  will  be  no 
reversal  of  the  search  direction  as  when  $ goes  past  the  solution. 
The  situation  will  be  similar  when  a near-tangency  condition 
exists  - i.e.,  when  the  curves  cross  twice,  but  the  crossing  points 
are  so  close  together  that  both  lie  between  successive  test  values 
of  <p . The  likelihood  of  this  happening  can  be  reduced  by  initially 
choosing  the  spacing  between  test  points  small  (i.e.,  choosing  the 
parameter  designated  as  N in  the  flow  chart  in  Figure  3.t-2  to  be 
large)  at  the  cost  of  lengthening  the  search  procedure.  In  any 
case,  the  algorithm  must  make  allowances  for  not  finding  a solution, 
if  only  to  accommodate  the  case  that,  due  to  errors  in  the  measure- 
ments, no  solution  compatible  with  the  measurements  exists. 
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FIGURE  3.^-2.  FLOW  CHART  FOR  ALGORITHM  TO  SOLVE  SINGLE-SITE 
GEOMETRY  BY  SEARCH.  TOA  T,  DAZ  a,  OAZ  b,  SLANT  RANGE  R, 
ALTITUDES  H AND  H ARE  GIVEN  MEASUREMENTS.  SEARCH  IS  SUB- 
ROUTINE (FIGURE  3?A-3) . N AND  M (FOR  SEARCH)  ARE  CHOSEN 
VALUES. 
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The  flow  charts  In  Figures  3.^-2  and  3.^-3  do  not  Indicate 
how  to  treat  cases  when  a solution  Is  not  found.  The  actual  solu- 
tion Is  likely  to  be  quite  sensitive  to  measurement  error  in  any 
case.  Thus,  unless  the  measured  slant  range  is  small,  it  may  not 
be  worth  finding  until  a new  set  of  measurements  is  obtained.  Two 
techniques  of  finding  the  solution,  wnich  may  be  used  in  combina- 
tion, are  to  decrease  the  size  of  the  increment  of  <t>  used  in  the 
search  and  to  change  slightly  the  measurement  of  R used  in  the 
calculation.  If  R is  changed  in  a direction  to  cause  the  passive 
radar  locus  curve  and  the  active  radar  locus  circle  to  intersect  if 
they  do  not  (and  to  intersect  at  two  more  widely  separated  points 
if  they  intersect  already),  then  the  search  algorithm  described 
will  be  more  likely  to  find  the  solution  of  the  modified  system. 

The  midpoint  between  the  solutions  using  the  modified  measure- 
ment set  should  be  approximately  the  same  as  the  midpoint  between 
the  solutions  based  on  the  original  measurements  (if  such  solutions 
exist)  and  should  be  a reasonable  approximation  to  the  actual 
target  position  if  no  solution  compatible  with  the  measurements 
exists  because  of  errors  in  the  measurements. 

These  extensions  to  the  algorithm  are  not  shown  in  the  flow 
charts  in  Figures  3.4-2  and  3- 4-3-  The  flow  charts  are  further 
simplified  in  that  they  omit  certain  details  that  should  be  in- 
corporated in  any  code  that  might  be  written  but  would  only 
obscure  the  basic  logic  if  included  here.  In  particular,  provi- 
sion should  be  made  to  prevent  division  by  zero  (e.g.,  if  a = 0) 
and  to  prevent  false  termination  of  the  search  because  the  test 
values  of  4>  appear  beyond  the  solution  region  ( | <t>  | > 180°  - |a|) 
due  to  numerical  roundoff  errors. 
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3.5  ACCURACY  OP  THE  SINGLE-SITE  SOLUTION 


A test  was  run  in  which  RADAR  was  located  20  nautical  miles 
due  south  of  OWN,  OWN  had  height  20,000  feet,  and  OTHER  flew  a 
level  semi-circle  of  radius  3 nautical  miles  centered  1000  feet 
above  OWN. 

The  bearing  of  OTHER  from  OWN  was  initially  taken  as  1 degree 
and  was  incremented  by  2 degrees  up  to  179  degrees.  The  computa- 
tion was  performed  10  times  at  each  step  with  random  normally 
distributed  input  errors  based  upon  sigmas  of  .15  microseconds  for 
the  two  TOA's,  100//12  feet  for  the  two  altitudes,  and  .2  degrees  for 
locating  the  center  mark  of  the  main  beam  passing  through  a point 
(either  North,  OWN,  or  OTHER). 

The  algebraic  solution  was  used.  Let  0 denote  the  bearing  of 
OTHER  from  OWN  relative  to  LOP.  Two  solutions  were  generally 
found  for  1°  < 9 < 65°  and  113°  < 6 < 179°.  When  1°  < 9 < 10°  and 
when  163°  < 0 < 179°,  sets  of  perturbed  inputs  that  produce  no 
solutions  are  not  uncommon.  It  is  also  not  uncommon,  especially 
in  the  second  of  these  regions,  for  the  perturbed  inputs  to  produce 
two  nearby  computed  values  of  <J>,  one  on  either  side  of  the  true 
value.  This  suggests  near  tangency  of  the  circle  determined  by 
active  interrogation  and  the  locus  curve  determined  from  the  passive 
measurements . 

The  error  sensitivity  (for  the  closest  computed  configurations, 
assuming  there  is  at  least  one)  varies  with  9.  In  the  best  range, 

9 not  very  far  from  90°,  the  standard  deviation  in  the  error  for 
9 is  a little  over  1/2  degree.  For  9 in  the  worst  ranges,  close 
to  0 or  180  degrees,  the  standard  deviation  is  about  2 degrees. 
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In  the  best  range,  the  standard  deviation  of  the  error  in  the 
computed  radar  distances  is  between  0.5  and  1 nautical  miles.  In 
the  bad  ranges  there  are  standard  deviations  of  over  5 nautical 
miles . 

There  is  large  inherent  error  sensitivity  when  the  absolute 
value  of  the  differential  azimuth  is  small  (i.e.,  when  OWN,  RADAR 
and  OTHER  are  nearly  colinear).  Given  such  a configuration,  it 
is  not  uncommon  for  the  algorithm  to  produce  no  solutions  from 
perturbed  inputs. 

The  algorithm,  as  coded,  writes  out  a default  solution  if 
either 

1)  no  solution  is  produced  by  the  general  method  or 

2)  the  differential  azimuth  is  so  small  that  error  sensitiv- 
ity makes  the  computed  solutions  suspect.  In  this  case 
the  computed  solutions  are  also  output. 

The  default  solution  involves  placing  OTHER  either  in 
between  or  opposite  RADAR  from  OWN.  The  choice  is  made  by  con- 
trasting the  TOA  as  measured  from  RADAR  with  OWN  TOA . (If  1/2 
OWN  TOA  is  less  than  TOA,  OTHER  is  opposite.) 

If  the  default  solution  is  used,  d is  not  computed. 


PASSIVE  MODE  RANGE-BEARING  CALCULATION: 

RADARS  WITH  AZIMUTH  REFERENCES 

il.l  OVERVIEW 

It  is  assumed  that  for  the  purpose  of  facilitating  the  imple- 
mentation of  BCAS,  modifications  will  be  made  to  either  all  or  a 
large  fraction  of  the  SSR  ground  sites  so  that  they  will  emit 
signals  allowing  ECAS  systems  to  determine  the  angular  orientation 
of  the  SSR  antennas  at  any  instant. 

This  will  allow  the  simplest  of  the  purely  passive  BCAS  opera- 
ting modes.  In  areas  where  there  is  coverage  by  two  suitably  placed 
radars  with  azimuth  references,  BCAS  can  determine  range  and 
tearing  to  a single  target. 

There  are  configurations  of  two  radars  and  a single  intruder 
for  which  the  solution  is  extremely  sensitive  to  perturbations 
of  the  input  variables,  i.e.,  the  measurements.  In  addition,  it 
will  be  seen  that  under  some  circumstances  measurements  from  two 
radars  will  give  multiple  solutions,  each  of  which  corresponds  to 
a distinct  radar-aircraft  configuration,  and  all  of  which  are 
completely  consistent  with  all  the  measurements  from  the  two  radars. 

In  either  of  these  troublesome  cases,  there  is  advantage  to 
using  ether  measurements  if  these  are  available.  One  may  then 
obtain  a more  accurate  target  position  determination  by  getting  a 
less  error-sensitive  solution  or  at  least  having  an  average  of 
several  error-sensitive  solutions.  When  the  measurements  from  two 
radars  result  in  multiple  solutions,  either  active  interrogation 
or  the  measurements  from  some  additional  radar  will  generally  be 


inconsistent  with  all  but  one  of  these  solutions,  which  can  then 
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be  taken  to  correspond  to  the  actual  configuration.  The  data 
management  schemes  appropriate  to  select  measurement  sets  and 
check  solutions  for  consistency  are  considered  to  be  beyond  the 
scope  of  this  report. 

The  particular  method  of  implementing  the  azimuth  reference 
signals  is  not  important  for  the  analysis  that  follows.  Similarly, 
the  question  of  whether  true  north  or  local  magnetic  north  should 
serve  as  the  reference  direction  is  not  considered.  In  the  simula- 
tions, it  is  assumed  that  the  reference  directions  for  the  two 
sites  are  the  same,  and  the  determination  of  own  azimuth  is  subject 
only  to  random  measurement  errors. 

Since  radar  altitudes  are  unknown,  the  radars  are  assumed  to 
be  at  sea  level.  The  nature  of  the  locus  curves  is  such  that  this 
assumption  should  have  no  significant  effect  on  the  results.  In 
general,  the  influence  of  small  altitude  variations  on  errors  in 
computed  positions  was  found  to  be  small. 

4.2  THE  SOLUTIONS  DESCRIBED  BY  THE  LOCUS  CURVES 

The  nature  of  the  solutions  to  the  problem  can  be  described 
in  terms  of  the  locus  curves.  In  addition,  an  enumeration  of  the 
possible  cases  also  suggests  that  the  solution(s)  can  be  found  by 
a search  algorithm.  The  algorithm  would  have  to  be  somewhat  more 
complex  than  that  for  the  single-site  geometry  described  in  Section 
3.4,  in  that  it  would  have  to  perform  search  in  two  dimensions, 
rather  than -one.  The  precise  choice  of  the  directions  for  search 
is  unimportant  as  long  as  the  whole  solution  region  can  always  be 
spanned . 
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The  principle  of  the  search  algorithm  would  be  similar  to  that 


in  Section  3.1*  - i.e.,  a locus  curve  lies  between  a pair  of  points 
such  that  T - Tc  changes  sign  in  going  from  one  to  the  other.  A 
solution  point  is  at  the  intersection  of  two  locus  curves,  or  the 
point  where  four  different  regions  (characterized  by  the  algebraic 
signs  of  - T^c  and  - T2c)  meet.  Thus  starting  from  one  of 
the  straight  lines  bounding  a solution  region,  the  search  algorithm 
would  find  a locus  curve  by  going  in  some  trial  directions  and 
then  find  the  intersection  by  searching  in  a direction  generally 
along  the  locus  line  already  found.  If  all  possible  cases  are 
properly  considered,  all  solution  points  can  be  found  to  within 
arbitrary  precision.  Difficulties  will  arise  if  near-tangencies 
occur  between  locus  lines.  The  solutions  in  such  cases  are  bound 
to  be  intrinsically  inaccurate  (measurement  error  sensitive),  but 
their  approximate  positions  can  be  determined,  as  in  the  single- 
site case,  by  modifying  the  measurement  parameters  (generally  the 
TOA's)  in  a direction  such  as  to  increase  the  intersection  region 
that  is  small. 

The  solution  regions  to  be  considered  are  classified  in 
Figure  4.2-1.  Knowledge  of  the  signs  of  the  differential  azimuths 
by  itself  restricts  the  target  position  to  one  of  the  four  regions 
indicated  in  the  figure.  The  regions  are  separated  from  each 
other  by  the  extended  lines  of  position  (LOP's)  between  OWN  and 
each  radar.  These  are  the  lines  passing  through  the  position  of 
OWN  whose  direction  is  determined  by  the  own  azimuth  (OAZ)  measure- 
ments . 
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FIGURE  it.  2-1.  DEFINITION  OF  SOLUTION  REGION  TYPES. 


The  regions  A,  B1  and  B2,  and  C are  qualitatively  different, 
and  the  nature  of  the  solution  in  each  may  be  different  when  the 
solution  is  sought  in  3 dimensions. 

Region  A is  bounded  by  the  extensions  of  the  LOP's  beyond  OWN 
away  from  the  radars.  It  is  characterized  by  having  no  radar  on 
its  boundaries.  Barring  exceptionally  unfavorable  situations 
(to  be  discussed  below),  the  locus  curves  for  solutions  in  this 
range  will  be  like  those  in  Figure  Jj.2-2.  There  will  be  one 
solution,  constrained  to  lie  in  the  sector  bounded  by  lines  in 
the  directions  6-^  + a^(a^<0) , and  6^  + a 2 (a2>0)  from  own.  The 
solution  point  can  be  determined  either  by  some  suitable  algorithm 
based  on  the  analytical  properties  of  the  functions  involved  or 


FIGURE  4.2-2.  USUAL  NATURE  OF  SOLUTION  IN  REGION  A. 

by  a search  algorithm.  Given  any  initial  point  in  the  region  of 
possible  solutions,  the  direction  of  search  for  an  improved  approx- 
imation can  be  determined  readily  from  the  sign  of  (T^-T^c)  and 

(T  -T  ) . 
v l2  2c 

Figure  4.2-3  shows  a potential  situation  when  there  might  be 
multiple  solutions  in  region  A.  This  might  occur  if  the  radars 
are  almost  directly  on  opposite  sides  of  OWN,  one  gives  rise  to 
locus  curves  with  two  branches,  and  the  other  has  a very  small 
differential  azimuth.  The  case  is  shown  (with  all  the  small  angles 
drawn  large,  for  clarity)  for  mathematical  completeness;  it  is 
highly  unlikely  in  practice.  Solution  1 is  the  usual  solution  in 
Region  A,  solution  2 the  alternative.  For  realistic  sets  of 
measurement  values  there  can  not  be  two  symmetrical  solutions  like 
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Solution  2.  The  proof  follows  from  arguments  like  those  used  to 
show  that  there  can  not  be  three  solutions  for  the  single-site 
geometry . 


FIGURE  2- 3.  EXCEPTIONAL  CASE  OF  MULTIPLE  SOLUTION  IN  REGION 
A (ANGLES  DRAWN  EXAGGERATED). 

Region  and  are  sectors  bounded  on  one  side  by  the  LOP 
from  OWN  to  a radar  and  beyond  and  on  the  other  by  the  extension 
of  LOP  beyond  OWN  away  from  the  radar.  They  are  characterized  by 
having  one  radar  on  their  boundaries.  Again  barring  exceptionally 
unfavorable  configurations,  the  locus  curves  determine  solutions 
according  to  the  scheme  of  Figure  4.2-1. 

It  is  seen  that  solutions  in  region  are  constrained  to 
lie  in  the  sector  bounded  by  lines  emanating  from  own  and  going  in 
the  directions  (B-^— it)  and  • Two  solutions  may  occur  in 
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then  either  the  "inner"  branch  of  the  radar  1 locus  curve  shown  in 
Figure  4.2-4  is  not  present,  or  it  does  not  intersect  the  locus 
curve  due  to  the  measurements  from  radar  2.  If  T,  > T,  every- 
where  along  the  LOP  from  own  toward  radar  1,  then  the  nature  of 
the  locus  curve  determined  by  the  measurements  from  radar  1 is  like 
that  shown  in  Figures  4.2-5  and  4.2-6.  The  question  of  the  exist- 
ence of  a second  solution  in  the  region  can  then  be  determined  by 
evaluating  T^c  and  T2c  along  the  other  boundary  of  the  region  in 
which  solutions  are  possible,  i.e.,  the  line  in  (usually)  the 
direction  from  OWN. 

Figure  4.2-7  shows  that  the  boundary  can  sometimes  be  the  line 
from  OWN  in  the  direction  8^  + “l*  If  a locus  curve  is  of  the  form 
in  Figure  2.2-2c  rather  than  of  the  form  2. 2- 2b  - i.e.,  not  going 
to  infinity,  the  nature  of  the  solutions  that  may  be  obtained  should 
generally  not  change . 

A special  case,  not  illustrated,  may  occur  if  both  radars  are 
in  essentially  the  same  direction  from  OWN  and  the  locus  curves 
have  "inner  branches"  that  intersect.  Then  there  will  in  theory 
be  three  solutions.  In  practice,  all  the  solutions  will  arise 
from  the  crossings  of  nearly  parallel  lines  and  will  be  extremely 
error  sensitive.  (This  is  intuitively  clear  - if  there  are  two 
radars  and  a target  nearly  colinear  with  OWN,  the  (small)  TOA  and 
DAZ  measurements  will  not  permit  good  determination  of  range  to 
the  target . ) 
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FIGURE  4.2-7.  POSSIBLE  UNUSUAL  BOUNDARIES 
OF  SOLUTION  REGION  B, . 


Region  C is  that  region  bounded  by  the  LOP's  - i.e.,  the 
sector  "between  the  radars"  when  viewed  from  OWN.  Given  that 
there  is  a solution  in  this  region,  there  may  be  one,  two,  or 
three,  depending  on  which  of  the  "inner  loops"  of  the  locus  curves 
exists  and  intersects  the  outer  loops,  as  indicated  in  Figure 
4.2-8  (a  theoretical  worst  case  is  five  solutions,  if  the  radar 
LOP's  are  so  close  that  the  "inner  loops"  also  cross  twice.  This 
is  of  no  practical  concern).  The  solution  indicated  as  Solution  1 
will  always  occur,  the  others  may  or  may  not. 

Figure  4.2-8  shows  another  possible  configuration  of  solutions 
in  Region  C which  may  occur  if  the  target  altitude  is  above  that 
of  OWN.  If  the  solution  exists  at  all  with  locus  curves  of  this 
nature,  then  either  the  pairs  of  solutions  1 and  2 or  the  pairs  1 
and  3 must  exist  together.  Barring  tangency  conditions,  there  are 
either  2 or  4 solutions.  Thus,  for  instance.  If  Solutions  1 and  2 
exist,  then  either  3 and  4 both  exist  or  neither  exists. 

Figure  4.2-10  shows  the  nature  of  the  solution  if 

T2  > 2 (H-Hq)  > Tr 

It  is  seen  that  generally  Solution  1 and  2 will  both  exist. 
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U.3  RANGE  AND  BEARING  COMPUTATION  ALGORITHM 


This  section  describes  the  computer  algorithm  developed  at 
TSC  and  Implemented  as  a FORTRAN  routine  to  run  on  the  DEC  PDP-10 
time-sharing  system  at  TSC.  The  algorithm  converges  to  a correct 
solution  except  in  some  cases  when  the  solution  is  extremely 
sensitive  to  perturbations  in  the  parameters  (i.e.,  measurement 
errors).  Since  the  algorithm  starts  by  finding  an  approximate 
solution  based  on  the  assumption  that  all  altitudes  are  zero, 
the  solution  generally  found  is  that  closest  to  the  planar  solu- 
tion. As  presently  implemented,  the  algorithm  does  not  attempt  to 
find  more  than  one  solution.  The  algorithm  Is  developed  in  the 
following : 

Assume  that  the  environment  of  OWN  contains  two  radars  at  sea 
level,  equipped  with  azimuth  reference  signals,  and  a target  air- 
craft . 

Let  r denote  the  horizontal  distance  from  OWN  to  the  target 
and  let  0 be  the  bearing  of  the  target  from  OWN,  relative  to 
North.  For  i = 1,2  let  d^  denote  the  horizontal  distance  from  OWN 
to  the  i-th  radar.  The  problem  is  to  determine  r,0  and  the  d^ 
from  the  input  BCAS  (given)  data. 

Let 

Hq  = height  of  OWN 
H = height  of  OTHER 

and  for  i * 1,2  let 

8^  = bearing  of  OWN  from  the  i-th  radar 

= (bearing  of  OTHER  from  the  i-th  radar)  - (bearing 
of  OWN  from  the  i-th  radar),  i.e.,  the  differential 
azimuth  of  the  target  measured  from  the  i-th  radar. 
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For  1 - 1,2  let  s^  be  the  horizontal  distance  from  radar  1 to  the 


target.  Let 

ri  r "1 


T,  = ^ s2  + H2  + ^ v2  + (H-Hq)2  - |/d^  + 


4.3.1 


— (where  c Is  the  speed  of  light)  Is  the  differential  time  of 
arrival  (TOA)  measured  from  the  1-th  radar.  The  given  data  con- 
sists of  H^j  H,  8^ , &2 » ^1*  a2 * ^ i > Tj • 

Consider  the  following  diagram,  where  the  points  labelled  OWN  and 
OTHER  are  really  the  horizontal  projections  of  OWN  and  OTHER. 


OTHER 


The  angle  of  the  triangle  whose  vortex  is  at  OWN  is  B.  - 9 and  the 
angle  whose  vortex  is  at  OTHER  is  it  - (6^-0+a  ). 

By  the  law  of  sines 


r sin( B^-9+ou  ) 
di  sin  a. 


4 . 3. 2 


s^  may  also  be  computed  in  terms  of  oujB^r  and  0 from  the 
geometry  of  the  same  triangle. 

The  expressions  for  d1  and  s^  may  be  substituted  in  4.3.1. 
The  results  is  of  the  form 

Tj_  = F(r,6,o1,B1,H,Ho)  i = 1,2  4.3.3 


78 


4.3*3  is  a system  of  two  equations  which  must  be  solved  for 


the  unknowns  r,0.  Then  4.3*2  computes  the  d^ . 

The  method  of  solution  involves  an  initialization  followed 
by  an  iteration. 

4.3.1  Initialization 

The  Initialization  assumes  that  the  radar  distances  d.^  are 
large  compared  to  r,  H and  Hq.  It  then  finds  a function  f(r,0,ai} 
61>H,Ho)  such  that  for  each  i,  f (r , 0 ,0^ , ,R,Ho)  will  approximate 
F (r ,0 ,a. , 6. ,H,H  ) and  solves  the  system 

llO 

= f(rJ0,01,H,Ho)  u 3 -i4 

Let  d(r,0,ct,B)  = r~Sl~r  ~ 6~9  ^ and  ror  each  1 let 

sin  ot 

d^  = d(r , 0 ,a^ , B^ ) • (Note  that  this  definition  is  consistent  with 
4.3-2) 

The  law  of  cosines  and  some  algebraic  manipulation  show  that 

r2  - 2rd . cos(0-P4)  + H2  - H2 
1 1 o 

+ r2  - 2djr  cosCb-b^  + H2  Wd2  + H2 


F(r,0,a1,B15H>Ho)  = 


r~2 2 4-3-5 

+ Y r + ( H-Hq ) d 

If  it  is  assumed  that  d^  >>  r,  H,  Hq>  then  the  denominator  of 
the  first  term  of  F is  close  to  2 d^,  and 

f (r,0,B,H,Ho)  = - r cosCB-B^  +Jr2  + (H-Hq)2  4.3.6 

is  a good  approximation  of  F. 

Let 

A = 1 - cosCBj-B-l)  4.3.7 

V = ((VT2)2  + (4-2A)TxT2  - 2A(H-Ho)2^  4 . 3 - B 


I 
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and 


2 (T. +T„ ) 2 2 

12--2—  (1J-2A)T]T2  + (A  -2A)(H-HQr 


it. 3-9 


The  system  of  equations  it.  3. 4 has  two  solutions  (r,0).  The 
two  values  of  r are  given  by  V+W. 

V+w  is  chosen  if  the  target  lies  in  the  smaller  wedge  determined 

by  Radar  - OWN  - Radar2>  Otherwise,  V-W  is  chosen. 

Once  r is  known,  it  is  not  difficult  to  solve  for  9. 

The  initialization  uses  the  (r,9)  computed  above  to  com- 
pute numbers  ^ in  such  a way  that  the  correct  solution  of  the 
system 

3^  = f(r,e,B1,H,Ho)  i = 1,2 

will  be  a better  initialization  than  the  one  previously  com- 
puted. The  solution  of  this  system  is  the  initialization. 

ft. 3. 2 Iteration 

Let  M(r,0)  denote  the  2x2  matrix 


4.3.10 


M(r,0)  = 

_3F 

3r 

(r,0,a1,B1,H,Ho) 

30  (r’e >a1>31,H,Ho) 

pz-i  I $-* 
ro  l<rt) 

J 

(r ,0 ,a2 , B2 ,h,HQ) 

||  (r,0,a2,62,H,Ho) 

If,  after  the  i—  iterative  step,  the  approximate  solution 
is  (r.^,0^)  then  ( r ^ + 1 > 9 ^ + i ) ls  defined  by  the  matrix  equation 

Tl"F(ri’9i’al’6l’H,Ho) 

T2-F(r^,0^,a2,B2,H,Ho) 


i + 1 


9i+l 


0 . 


+ M ( r , 0 ) 


-1 


4.3.11 


The  iteration  stops  if 

1)  i = 20 

or 

2)  ri+i-rj[  and  9i+l“®i  are  sufficiently  small. 


R0 


The  iteration  is  along  a vector  field.  If  (r,0)  is  the 
solution,  then  (r,9)  is  a sink  of  this  field,  but  It  is  not 
the  only  sink.  In  particular,  (-r,9)  (which  is  geometrically 
meaningless)  and  (r,7r+0)  are  both  sinks.  If  either  of  these 
two  additional  sinks  is  chosen  by  the  iteration,  the  computed 
position  of  the  other  will  be  the  reflection  about  OWN  of  the 
true  position.  Fortunately,  it  is  not  difficult  to  tell,  from 
the  Input  data,  if  this  reflection  has  been  computed.  There- 
fore, the  algorithm  can  correct  the  "mistake"  of  finding  the 
wrong  sink. 

b.h  NATURE  OF  THE  SOLUTION 

The  algorithm  described  in  the  previous  section  converges 
to  a correct  solution  except  when  the  aircraft  and  radars  are 
in  certain  configurations  in  which  the  solution  is  extremely 
sensitive  to  small  variations  in  the  measurements . The  iteration 
then  may  not  converge. 

The  term  "correct  solution"  designates  a solution  which 
describes  a configuration  for  which  the  noise-free  TOA,  DAZ 
and  OAZ  measurements  would  be  precisely  those  from  which  the 
solution  was  obtained.  Generally,  only  one  such  solution  exists, 
but  under  some  circumstances  there  may  be  several  solutions. 

Such  a case  is  illustrated  by  the  third  simulation  run  described 
in  Section  5 and  is  discussed  below. 

It  is  important  to  state  again  that  the  false  solution  is 
not  a construct  due  to  some  deficiency  in  the  algorithm  or  to 
excessive  measurement  errors. 


It  is  instead  one  of  two  configurations  of  radars  and 
aircraft  which  result  in  identical  sets  of  measurements  (of  i 

TOA,  DAZ,  OAZ  and  altitudes)  being  made  by  the  BCAS.  When  two 
distinct  configurations  are  possible  sources  of  the  measurements 
obtained,  the  nature  of  the  algorithm  determines  which  it  will 
find,  but  it  is  essentially  a matter  of  chance  whether  this  is 
in  fact  the  actual  configuration. 

Though  calculations  made  at  only  one  instant  and  based  only 
on  the  nominal  inputs  cannot  determine  whether  a particular  solu- 
tion that  has  been  obtained  corresponds  to  the  actual  configu- 
ration, in  practice  there  may  often  be  additional  data  available 
that  allow  the  decision  to  be  made.  It  should  be  possible  to 
extend  the  algorithm  to  determine  the  alternative  solutions  when 
these  exist. 

A number  of  checks  can  be  performed  on  the  solution:  If 
the  range  to  one  or  both  of  the  radars  has  been  calculated  in 
some  other  context,  a rough  consistency  check  can  be  performed. 

Active  interrogations  can  be  performed  to  get  an  independent 
estimate  of  the  range  to  the  target  as  a check  on  the  value  com- 
puted from  the  passive  measurements.  Finally,  if  a time  sequence 
of  successive  configurations  is  computed,  though  the  improper 
solution  may  be  found  at  each  time,  these  configurations  will 
not  be  consistent  in  time.  Thus,  in  the  case  of  the  example  in 
the  simulated  sequence,  it  may  be  observed  that  the  calculated 
distance  between  the  radars  is  about  38  miles  when  the  improper 
solution  is  first  found  and  continuously  changes  to  about  56 


miles  some  46  seconds  later.  It  must  be  noted,  however,  that 
the  calculated  distance  to  the  radars  is  quite  sensitive  to 
measurement  errors  and  shows  large  random  variation  between 
successive  instants  of  calculation.  Thus  inconsistencies  in 
calculated  distances  to  radars  can  be  taken  to  be  reliable  in- 
dicators of  an  improper  solution  only  if  they  are  very  large  or 
sustained  over  a period  of  time. 

A more  satisfactory  answer  to  the  problem  of  multiple  solu- 
tions would  be  an  algorithm  that  determined  all  solutions  wher- 
ever several  exist.  Further  software  in  the  system  could  then 
resolve  the  ambiguity,  either  by  calling  for  active  range  mea- 
surements or  by  evaluating  the  track  data.  It  may  be  anticipat- 
ed that  when  several  distinct  solutions  exist  which  are  quite 
different,  it  should  be  easy  to  determine  the  actual  configura- 
tions. When  the  solutions  are  distinct  but  relatively  similar, 
the  resolution  of  the  ambiguity  may  be  more  difficult. 

In  addition  to  configurations  for  which  the  BCAS  measure- 
ments are  such  that  several  distinct  consistent  solutions  may 
be  calculated,  there  exist  configurations  in  which  relatively 
large  displacements  of  the  targets  or  radars  lead  to  small 
changes  in  the  measurements  made  by  BCAS.  Inversely,  small 
changes  in  the  measurements  (such  as  those  due  to  quantization 
and  noise)  lead  to  large  changes  in  the  configuration  corres- 
ponding to  the  measurement  set.  This  means  that  the  values  of 
range  and  bearing  calculated  from  noisy  measurements  will  show 
large  deviations  from  their  actual  values  as  a result  of  the 
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random  errors. 


Two  problems  then  arise  in  these  regions:  First,  the 
values  of  range  and  bearing  that  may  be  found  when  the  equations 
describing  the  measurements  are  solved  may  be  quite  far  from 
the  actual  values.  Second,  the  nature  of  the  problem  is  such 
that  the  iterative  algorithm  described  above  for  solving  the 
system  may  fail  to  converge. 

The  "bad"  configurations  generally  are  the  following: 

1.  Both  radars  in  the  same  direction  from  the  BCAS  air- 
craft (OAZn-OAZ2  less  than  'vlO0). 

2.  The  intruder  aircraft  between  the  BCAS  and  one  of  the 

ground  radars  ( | —9  | less  than  'v  15°). 

3.  The  intruder  aircraft  in  the  direction  opposite  the 
radar  from  BCAS  (( (5^.-0-  180°  j less  than  ^ 3°)> 

The  precise  extent  of  the  bad  ranges  depends  also  on  the 
other  aspects  of  the  configuration — altitudes,  distances  to  the 
radars,  distance  of  the  target,  and  the  magnitude  of  the  bearing 
angle  that  is  not  the  primary  contributor  to  the  problem.  In 
general,  the  bad  ranges  tend  to  be  wider  when  the  aircraft  are 
at  higher  altitudes.  In  the  case  of  configurations  2 and  3S 
when  the  intruder  is  colinear  with  the  BCAS  and  the  radar,  the 
bad  range  widens  and  splits  when  the  aircraft  are  at  high  alti- 
tudes. The  errors  in  the  computed  positions  are  large  when  the 
intruder  aircraft  is  on  either  side  of  the  radar-BCAS  line,  but 
moderate  when  directly  on  the  line. 


All  of  the  enumerated  effects  can  be  explained  qualitatively 
by  the  behavior  of  the  locus  curves  shown  In  the  Figures  In 
Section  2. 
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5,  PASSIVE  MODE  BCAS  WITHOUT  AZIMUTH  REFERENCE  SIGNALS 

5 . 1 OVERVIEW 

A purely  passive  BCAS  system  without  azimuth  reference 
signals  can  operate  when  there  are  two  radars  and  two  targets. 

Since  all  the  measurements  available  to  the  system  are 
invariant  under  rotation  of  the  whole  configuration,  the  bearing 
of  one  radar  may  be  specified  arbitrarily.  The  three  relative 
bearing  angles  in  the  configuration  are  defined  in  terms  of  the 
arbitrary  reference  direction. 

The  quantities  measured  are  four  TOA's,  four  DAZ's,  and 
three  altitudes,  or  eleven  measurements  in  all. 

The  interrelationship  between  measured  and  derived  parameters 
is  shown  in  Figure  5.1-1. 

Solutions  to  be  found  are  for  the  parameters  describing  the 
configuration.  These  are  two  distances  to  the  ground  interroga- 
tors from  OWN,  two  distances  to  both  targets,  three  bearing  angles 
and  three  aircraft  altitudes,  or  ten  parameters  in  all.  Thus, 
the  problem  is  overdetermined. 

The  method  of  solution  is  complicated  and  consists  of  a 
number  of  steps.  At  the  end,  a solution  is  sought  which  in  a 
useful  sense  is  a best  fit  to  the  data.  The  solution  proceeds  in 
three  separate  essentially  independent  steps  - coarse  initializa- 
tion, iterative  refinement,  and  least-squared-error  fitting.  In 
outline  the  steps  are  the  following. 

Step  1:  The  coarse  initialization  is  based  on  the  simplify- 
ing assumptions  that  the  BCAS-to-target  distances  are  much 
smaller  than  the  BCAS-to-radar  distances  and  that  all  aircraft 


86 


(b) 


FIGURE  5.1-1.  DEFINITION  OF  VARIABLES  IN  TWO-RADAR,  TWO 
TARGET  CONFIGURATIONS.  (a)  SHOWS  HORIZONTAL  DISTANCES, 
(b)  DEFINES  ANGLES  ABOUT  OWN.  ALL  ANGLES  SHOWN  ARE 
POSITIVE  IN  SIGN. 
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altitudes  can  be  neglected.  These  simplifications  allow  the 
inherently  non-linear  set  of  eight  equations 


Tij  = Sij  + RJ  “ Di 


5.1.1 


2 2 

3 + d - r 

cosoi. , = — ^ 

lj  2s44  d4 


5.1.2 


ij  ui 


(1=1 > 2 ; J=l,2) 


to  be  reduced  by  various  algebraic  manipulations  and  successive 
elimination  of  variables  to  a pair  of  simultaneous  linear  equations 
that  is  solved  for  the  target  distances  r^  and  r^.  The  original 
set  of  equations  for  this  simplified  case  is  still  overdetermined 
(assuming  the  altitude  to  be  0,  there  are  8 equations  correspond- 
ing to  the  measurements  to  be  solved  for  seven  parameters  defining 
the  configuration).  Therefore  inconsistent  sets  of  values  can  be 
obtained  for  the  other  variables  defining  the  configuration, 
depending,  on  the  order  in  which  the  variables  already  solved  for 
are  substituted  back  into  the  equations  to  evaluate  the  others. 

The  choice  made  is  to  determine  each  radar-to-BCAS  distance  on  the 
basis  of  the  greater  DAZ  (so  as  to  minimize  the  percentage  error 
due  to  measurement  errors).  The  other  quantities  are  computed  in 
an  essentially  random  order,  with  no  attempt  made  to  minimize  the 
resultant  inconsistencies. 

Step  2:  The  coarse  initialization  obtained  in  Step  1 serves 
as  the  initial  configuration  to  be  improved  iteratively  in  Step  2. 
The  main  goal  of  the  process  is  to  take  into  proper  account  the 
aircraft  altitudes,  neglected  during  Step  1.  Further,  explicit 
note  is  taken  of  the  fact  that  the  system  of  equations  is  over- 


determined . 


The  calculations  are  performed  iteratively.  Thus  distances 
and  angles  are  computed  for  the  projection  of  the  3-dimensional 
configuration  of  the  aircraft  on  the  horizontal  Diane.  The  effect 
of  the  aircraft  altitudes  is  taken  into  account  by  modifying  the 
measured  values  of  TOA  to  compensate  for  the  altitude. 

One  may  note  that  the  TOA 

T = /s2  + H2  + 2 + (H-Hq)2 

- Vd2  + H2  5.1.3 

= s + + r + eR  - d - eD 

where 

cs  * h2  - s 5.1.1* 

eR  =|/r2  + (H-Hq)2  - r 5.1-5 

eD  + Ho  ” d 5.1.6 


H - altitude  of  target 

H - altitude  of  BOAS 
o 

One  can  therefore  write  a simplified,  apparently  linear  TOA 
equation 


where 


t = s + r - d 


5.1.7 


t = T - es  - eR  + eD  ' 


5.1.8 


At  each  step  of  the  iteration,  the  value  of  t is  recomputed  on  the 
basis  of  the  best  current  estimates  of  s,  r and  d until  convergence 
is  achieved  - i.e.,  until  the  values  do  not  significantly  change 
from  step  to  step.  (In  certain  "bad"  configurations  the  process 
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does  not  converge  and  it  is  therefore  always  terminated  after  some 
fixed  number  of  iterations.) 

A set  of  three  independent  equations  in  the  two  unknowns  d^ 
and  d^  is  developed.  This  overdetermined  system  is  solved  at  each 
iterative  step  in  the  following  way: 

Three  equations  developed  on  the  basis  of  geometric  arguments 
require  that: 


F1(d1,d2)  = 0 

F2(d1,d2)  = 0 5.1.9 

P3(d1,d2)  = 0 


where  F^,  F2,  and  F^  are  functions  only  of  the  BCAS-to-radar 
distances  (d^  and  d2),  the  measured  differential  azimuths,  and 
(the  best  current  estimates  of)  the  adjusted  TOA's  (Equation  5.1.3)' 
These  equations  are  mutually  inconsistent  - i.e.,  for  any 
pair  of  values  (d^,d2),  not  all  three  functions  F^,  F2,  F^  will  be 
identically  zero,  but  rather  they  will  have  values 


such  that 


T?  ( h d ) = e 
1 1 VU1 >u2 ; 1 

F2(d1,d2)  = e2 

F3(d1,d2)  = e^ 

* (A  ♦ A * 4) 


5.1.10 


5.1.11 


At  each  step,  the  iterative  algorithm  determines  changes  (Ad^, 

Ad2)  to  d^  and  d2  such  as  to  reduce  E,  the  measure  of  the  incon- 
sistency of  the  equations. 

The  changes  (Ad^,Ad2)  are  computed  as  follows:  The  equations 
are  made  into  linear  equations  in  Ad^  and  Ad2> 
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(d^  + Ad^,  d^  + Ad2) 


5-1.12 


3F  3K, 

= F1  (d1,d2)  + ^d7  Adj  + jj-  A' 


= e1  (Ad1,  Ad2) 


One  seeks  the  values  of  Ad,  and  Ad_  which  minimize 


E (Ad^,  Ad2)  - Z (Ad.^  Ad2)‘ 


5.1.13 


The  minimum  occurs  when 


(j  = 1,2) 


5.1.1^ 


This  is  a set  of  two  linear  equations  in  the  two  unknown  Ad^  and 


Ad2-  Its  solutions  are  used  to  improve  the  current  values  of  d.^ 


and  d2,  to  compute  the  other  parameters  that  determine  the  con- 


figuration from  these,  and  to  update  the  values  of  the  adjusted 


TOA's.  The  iterative  step  is  then  repeated  until  convergence  is 


obtained  (or  failure  to  converge  is  evident). 


When  this  process  is  used  to  compute  the  configuration  from 


simulated  noise-free  measurements,  perfect  results  are  obtained 


(where  the  process  converges).  When  noisy  measurements  are  used 


(i.e.,  in  the  practical  case),  then  reasonably  good  fits  to  the 


actual  configuration  are  obtained.  However,  these  are  not  the 


best  fits  to  the  data.  Furthermore,  as  in  the  case  of  Step  1,  the 


parameters  defining  the  configuration  are  determined  by  first 


finding  one  pair  of  them  - here  d^  and  d?  - and  then  determining 


the  rest  successively  by  substituting  the  values  of  the  parameters 


already  solved  for  into  expressions  involving  the  others.  Since 


the  overall  set  of  equations  is  overdetermined,  the  values  obtained 


■^X’: 


L 


will  not  in  general  be  consistent.  No  attempt  is  made  in  this 
step  to  revolve  these  inconsistencies.  The  theoretically  optimum 
solution  is  obtained  by  Step  3,  Tor  which  the  results  obtained 
here  serve  as  initial  values. 

Step  ?.  The  final  step  is  again  an  iterative  squared-error 
minimization  process.  The  eleven  quantities  measured  by  BCAS  are 
expressed  as  functions  of  the  ten  various  coordinates  defining  the 
radar-aircraft  configuration. 

(X)  5.1.15 

where  y^  is  the  H-th  measurement  and  X = (x^  ...  x]_q)  is  the  vector 
of  coordinate  values  defining  the  radar  aircraft  configurations. 

The  components  of  X are  the  two  BCAS-to-radar  distances,  the  two 
BCAS-to-aircraft  distances,  the  three  relative  angles  to  aircraft 
and  radar  from  the  BCAS,  and  the  three  aircraft  heights.  The  y^ 
are  the  following:  four  TOA ' s , three  reported  altitudes,  and  for 
each  radar  the  sum  and  difference  of  the  differential  azimuths  of 
the  two  target  aircraft  relative  to  that  radar.  (The  sum  and  dif- 
ferences of  the  DAZ's  are  used,  rather  than  the  DAZ's  themselves, 
because  there  is  correlation  between  the  measurement  noise  com- 
ponents of  the  DAZ's,  but  not  between  the  noise  components  of  their 
sums  and  differences.) 

The  actual  measurements  m^  are  noise  corrupted,  so  that  there 
will  be  a random  discrepancy  e^  between  the  predicted  value  y^  of 
a given  measurement  when  the  configuration  is  described  by  a given 
set  of  parameters  X and  the  actual  measurement  m^. 

ez  = m^  - yz(X)  5-1.16 
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This  discrepancy  e is  ascribed  to  measurement  error.  By  what  is 
known  as  the  principle  of  least  squares,  the  assumed  configuration 
test  fits  the  measurement  data  when 


i 4 

l O £ 


5.1.17 


= I 

Si 


y,«))2 


is  minimized  (where  a^'s  are  the  variances  of  the  independent 
errors  in  the  measurements). 

The  X minimizing  E is  found  iteratively.  The  set  of  equations 
for  the  errors  are  first  made  into  a set  of  linear  equations  in 
terms  of  AX,  incremental  changes  about  the  true  minimum  configura- 
tion. The  set  of  equations  is  of  the  form 


Fg,  (Xopt^  _ y 

aH  k 9xk 


5.1.18 


The  partial  derivatives  are  evaluated  at  the  current  best  approxi- 
mation of  the  true  configuration.  Temporarily  holding  these 
partial  derivatives  fixed,  standard  multivariate  regression  tech- 
niques are  used  to  find  the  AX  that  minimizes  E.  The  set  of 
coordinates  X is  then  corrected  by  adding  the  computed  AX  and  the 
process  is  repeated  until  it  converges  - i.e.,  until  successive 
AX's  become  sufficiently  small. 

5.2  INITIALIZATION 

The  case  considered  is  illustrated  in  Figure  5.1-1.  The 
notation  throughout  is  the  following. 


r = the  horizontal  distance  from  OWN  to  target  j 
J 

= the  horizontal  distance  from  OWN  to  radar  i 

s^j  = the  horizontal  distance  from  radar  i to  target  j 

Rj  , Dj  and  S^j  are  the  corresponding  slant  distances 

a^j  is  the  differential  azimuth  of  target  j measured  from 
radar  i 

<t> ^ j is  the  angle  between  the  i-th  radar  and  the  j-th 
target,  viewed  from  OWN 
H is  the  height  of  OWN 


is  the  height  of  the  j-th  target 


For  i = 1,2,  j = 1,2 


T =J s?.  + U2.  +*/r?  + (F.-H  )2  ~%/a2  + H2 

ijfij  j y j j o y i o 


5.2.1 


= (bearing  of  target  j from  radar  i)  - (bearing  of  own  from 
radar  i.  From  the  law  of  consines 


rj  = sij  + di  " 2’Vsij  cosaij 


5.2.2 


Equation  5.2.1  is  the  (i,j)  TOA  equation  and  Equation  5.2.2  is  the 

(i,j)  differential  azimuth  equation. 

Assume  all  heights  are  zero  and  that  the  distances  r.  (BOAS 

J 

to  target)  are  small  compared  to  the  distances  d^  (BOAS  to  radar). 


By  the  law  of  cosines 


2 2 2 

s = r.  + d.  - 2r,d.  cos0, , 

ij  j i j i ij 


Hence,  by  the  TOA  equation 


rj  + di  - 2rjdi  coseij  + rj  - di 


5.2.3 


5-2.4 


1 
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a 


Vdl  * cos9lJ  * (rj/dj)'  - 1 ♦ 

By  the  approximation 


(1+x)  = 1 + nx  ( | x | <<1) 

which  will  be  used  repeatedly  in  what  follows 


Tij/di  = - rj/di  cos<^ij  + (r^/d^K  + rj//di 


r^/d^^  (l-cos^) 


Multiplying  by  d. 


Tij  = rj  (1-COS<frij) 


According  to  the  differential  azimuth  equation 


ri2  = sij  + di  - 2'sij*d1-cosctij 


Therefore 


cosa 


2 . ,2  2 
si.i  * di  - r.i 


U 2sij'di 


Substitution  of  5.2.1  in  5.2.10  and  simplification  yields 


cosa. . 


di  “ r1  cos<t>1.1 


iJ  l/Ti  * di  ■ co 


s4> 


iJ 


Let  t^j  = r^/d^ 


Then  5-2.11  becomes 


COSlj) 


1-t1.1  COS*li 


ij 


l"2tij  cos*ij  + fcij 


1“ ^ j COS^^j 


/ 2 5 2 

w(l-t,,  cosd)..)  + t..  sin  4> . . 


.2.5 

.2.6 

.2.7 

.2.8 

.2.9 

.2.10 

.2.11 


5.2.12 


and,  since  1-t^,  cos<J>1j  is  positive 


r 


1 


5.2.21 


? A 


ii 


Ti J r.i  - Tu 


Hence 


' rj/di 


ti2/,tll  " r2/rl 


5.2.22 


Divide  the  (1,2)  equation  of  5.2.21  by  the  (1,1)  equation  and 
use  5.2.20.  The  result  is 


2 r - T 
1 il 


<?  r - T 
^ r2  12 


5.2.23 


When  5.2.23  is  squared,  the  result  is  a linear  equation  in 
the  variables  r^ , r^.  There  are  two  such  equations  corresponding 
to  i = l , 2 . These  may  be  solved  simultaneously  for  r^ , r^ . 

Given  i,  choose  j such  that  Ajj  is  maximal  and  compute 

T. 


di  = ri/1:iJ 


1LL 


2 A 


U 


2 r,  - T 


ij 


5-2.2H 


s^,  may  now  be  computed  by 

Sij  = TU 


rj  + di 


5.2.25 


5.3  ITERATIVE  IMPROVEMENT 

This  section  describes  an  algorithm  that  takes  the  output  of 
the  initialization  algorithm  of  the  previous  section  and  produces 
a result  that  is  correct  (except  for  the  bad  ranges  described  in 
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(5.6)  if  the  measured  quantities  (TOA's,  differential  azimuths  and 
heights)  are  correct. 

An  observation  that  is  central  to  the  algorithm  is  that  the 
TO A equations 

Tu  -fAj + - J«l  * »o  5-3-1 

simplify  considerably  to 


T 


ij 


ij 


+ r. 

1 


d. 

i 


5.3.2 


if  the  heights  are  assumed  to  be  0. 

If  approximate  values  are  known  for  s^j , d^  and  r . the  error 
in  5.3.2  can  be  approximated  and  absorbed  into  the  left-hand  side 
of  the  equation. 


Let 


5.3.3 


Then  from  5-3.1  and  5-3-3 

' su  * q - ai  5.3.1 

Given  old  values  for  the  variables  d^,  r ^ , s^j  the  algorithm  uses 
5.3.3  to  approximate  the  U^.  and  then  uses  the  simplified  TOA 
equations  5.3.4. 

There  is  another  important  observation  to  be  made.  Super- 
ficially, it  seems  that  the  system  to  be  solved  consists  of  eight 
equations  (the  four  TOA  and  four  differential  azimuth  equations) 
in  the  eight  unknowns  d^  r . , s^. . There  is,  however,  an  additional 
constraint  that  must  be  placed  on  the  d^,  r.,  s_  so  as  to  ensure 
that  this  variable  set  actually  corresponds  to  a geometric 
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configuration.  One  way  of  seeing  this  is  to  count  dimensions. 

The  input  measurement  space  is  11-dimensional  (four  TOA's,  four 
differential  azimuths  and  three  heights).  The  configuration  space 
locates  four  points  (the  radars  and  targets)  on  the  ground,  eight 
dimensions,  and  includes  three  heights  for  a total  of  11,  but  two 
configurations  are  'the  same'  if  they  differ  only  by  a rotation 
about  own.  This  takes  away  one  degree  of  freedom  from  the  config- 
uration space.  Thus,  the  set  of  measurement  vectors  that  corre- 
spond to  geometric  configurations  is  a 10-dimensional  subset  of  an 
11-dimensional  set.  If  a measurement  vector  does  not  lie  in  this 
subset,  there  is  no  configuration  that  corresponds  to  that  vector. 
For  'most'  measurement  vectors,  then,  it  is  impossible  to  solve 
for  a configuration  and  reality  is  best  represented  by  an  over- 
determined system. 

There  are  two  obvious  ways  to  modify  the  original  system  of 
eight  equations  in  eight  unknowns; 

1.  Change  the  variable  set  describing  the  configuration  from 
{dj,  r j , s^j,  Hj } to  a set  of  ten  variables  in  such  a 

way  that  each  set  actually  corresponds  to  a configuration. 

2.  Add  an  equation  to  the  original  system  of  eight  equations 
in  such  a way  that  if  the  equation  is  satisfied  then  the 
variable  set  corresponds  to  a configuration. 

The  second  of  these  options  will  be  chosen. 


Let  U^j  be  defined  as  in  5.3.3  so 

Uu  * su  + rj  ■ 


su  * 


rj  * di 


5-3.5 


Substitution  of  5.3.5  into  the  ( i , j ) differential  azimuth 


5.3.9  represents  a system  of  two  equations  (j=l,2)  in  the  two 


unknowns  d^,  d0.  Once  these  two  variables  are  known,  it  is  easy 


to  find  the  others. 

It  was  mentioned  earlier  that  the  computation  of  the  variables 
dj,  r.,  s^j  involves  solving  an  overdetermined  system  of  equations. 
Therefore,  another  equation  must  be  found. 

<J>  . will  be  taken  to  be  a positive  angle  of  less  than  180 
degrees  if  and  only  if  is  positive.  With  this  convention,  the 
angle  with  vertex  0 that  is  spanned  by  d^  and  d2  is  given  by 

2J  ®1J* 

Thus 


0 = sin(4>22-$.l2)  - sin($21-$11) 


5.3.10 

5.3*10  can  be  written  in  terms  of  the  sines  and  cosines  of  the  $. 


1J 


and  will,  then,  yield  a third  question  in  d^,  d2  if  the  sin  , 
cos$^j  are  expressed  in  terms  of  d^,  d2> 

By  the  law  of  cosines 


cos$ 


U 


2 

. !u_ 


2 2 
di  'll 


2 di  r3 


Also 


and  5-3.7 


sij  = Uij  + di  - rj 


_ Al.l‘di  + A1.1U1.1  di  + 2 Ui,1 
J U„  + A,  , d. 


5.3.11 


5.3.12 


5.3.13 


U 


ij 


Substitution  of  5. 3-12  and  5. 3-13  in  5.3*11  gives  cos#^  as 
a function  of  d^,  d,,. 


5. 3. 1U 


Similarly,  by  the  law  of  sines, 


si  1 

sin$,  . = — >L  sina.  . 
1J  rs  ij 


and  substitution  of  5. 3-12  and  5-3.13  in  5-3-14  gives  sin^j  as  a 
function  of  d^,  d^- 

When  5-3-10  is  expressed  in  terms  of  sin<j>^  , cos<t>^  the 
result  is  a third  equation  in  d^,  d^ . 

The  iteration  may  now  be  described. 

An  iterative  step  starts  with  approximate  values  of  the  d^, 
r j , s^j  . Next  the  U„  are  computed  by  5-3-3- 
Let 

0 = F1(d1,d2) 

0 = F2(d1,d2) 

be  the  two  equations  of  5-3-9  and  let 
0 = F3(d1,d2) 

be  the  equation  derived  from  5-3-10 

The  iteration  proceeds  by  finding  an  approximate  solution  to 
the  system 

0 = F1(d1,d2)  5-3-15 

0 = F2(d1,d2) 

0 = F3(d1,d2) 

then  computing  the  rj  and  s^.  in  terms  of  the  d^^  and  then  going  to 

the  next  iterative  step  unless  the  changes  in  the  variables  have 
been  small  enough. 
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It  only  remains  to  describe  how  the  approximate  solution  to 


5.3.15  Is  computed. 

Suppose  d^,  d^  denote  the  approximations  obtained  prior  to 
the  current  iterative  step.  The  object  is  to  find  Ad^,  Ad,,  in 
such  a way  that  d^  + Ad^,  d,,  + Ad2  are  better  approximations. 

The  best  linear  approximation  to  5.3.15  Is  given  by  the 
matrix  equation 

1 ^pl  J}T?1 

- FJ'(d1,d2)  = (dltd2)  Adx  + (d1,d2)  Ad2  5.3.16 

- p2fdl.d2>  ■ (dl-d2>  4dl  + $£;  <dl-d2>  4d2 

- f,3<d1-d2>  - <dl-d2>  4dl  * % (dl>d2)  4d2 


th 


The  approximate  solution  to  5. 3. 16  that  is  found  is  the 
weighted  least  square  solution  where  the  weights  are  chosen  as 
follows : 

First,  each  equation  is  weighted  in  such  a way  that  the  i 
row  of  the  coefficient  matrix  of  5.3-16  will  have  norm  1. 

For  j = 1,  2,  let  r^  = min  (U^j ,U2j ) . 
t h 

The  j equation  of  the  normalized  system  is  then  weighted  by 

r . 


+ r; 


, for  j = 1,  2 


No  claims  are  made  for  the  optimality  of  this  set  of  weights. 
They  were  chosen  by  trial  and  error  and  probably  can  be  improved 
on. 
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5.4  LEAST  SQUARES  FITTING 

5.4.1  Derivation  of  the  Normal  Equations 

Assume  that  there  are  defined  M functions  y , which  are  all 
linear  combinations  of  N independent  variables  xn,  where  M > N. 

N 

y = V 6 x 5.4. 

3 m mn  n 


The  8 are  taken  to  be  known  constants, 
mn 

as  noise-corrupted  measurements. 


? 


The  ym  are  given 

5.4.2 


where  nm  is  the  noise  component  in  the  m-th  measurement,  assumed 
to  be  random  with  zero  expected  value.  The  problem  is  to  deter- 
mine the  values  of  the  *n's  that  best  explain  the  measurements 
obtained.  (This  problem  is  somewhat  different  in  principle  from 
a more  familiar  problem,  in  which  both  the  ym  and  a set  of 
variables  xmn  are  measured,  and  the  task  is  to  determine  a set 
of  & defining  the  functional  relationship  of  the  measurements. 
However,  the  equations  to  be  solved  will  have  the  same  form, 
with  the  3's  and  x's  in  interchanged  roles.) 

A measure  of  how  well  the  x^s  explain  any  given  measurement 

y ' is  the  residual 
m 

M 

r = ym  - X &mn  5.4.3 

m m , mn  n 
n=l 

The  principle  of  least  squares  states  that  the  best  overall  explana- 
tion y'  is  that  set  of  x 's  for  which 
m n 
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1 


R(x,  x0...x  ) = l,  r 

1,  2 n m=l  m 

is  a minimum.  The  conditions  satisfied  at  the  minimum  are  that 


= 0 (n  = 1,2. . .N) 


5 - ^ ^ 


M 3 r 

E - 

m=l  3 x 


5-4.5 


I,  2 r 

iri”  1 m « 

3xn 


5-^-6 


M N 

_3R_  , 2 2(y^  - Z 3 £ x )(-Bmn) 
3x  m=l  m S.=  l 


5.4.7 


M N M 

= _2rr^l6nnvm  +2^1  xl 


Hence,  there  are  N equations 


M N M 

m=l  ^mn  ym  ~ £-l^m=l  ^mJl^mn^  xi 


5.4.8 


to  be  solved  for  the  values  of  x^.  These  are  N linear  equations 
in  N unknowns,  and  the  solution  is  straightforward.  Using 
matrix  notation,  if  corresponding  to  equation  5.4.1 


Y * BX, 


5.4.9 


then  the  solution  to  5.4.8  is 


X = (B  B1)  B1  Y . 


5.4.10 
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It  is  a known  result  that  when  the  errors  n are  normally 

m 

distributed.  Independent,  and  with  zero  means  and  equal  vari- 
ances, the  resultant  solution  is  also  a maximum  likelihood  esti- 
mate of  X. 


5.4.2  Linear  Approximation  of  the  BCAS  Equations 

The  equations  relating  the  positions  of  the  aircraft  and 
radars  to  the  quantities  measured  by  BCAS  are: 

The  TOA  equation 


* 'W 


5.4.11 


The  DAZ  equation 


°i,J=  - cos 


-1 


rl 


‘disU 


5.4.12 


where  the  sign  must  be  chosen  properly. 

The  subscript  i designates  the  radar  and  j the  target. 

Tj  is  the  horizontal  distance  of  OWN  to  target;  d^  is  the 
horizontal  distance  of  OWN  to  radar;  H and  H,  are  OWN  and 
target  altitudes;  s^j  is  the  distance  of  target  to  radar. 

si,j  = (di2  + rj2  ~ 2dirjcos0i,J}  5.4.13 

where  ^ is  the  angle  at  the  BCAS  between  the  direction  toward 
the  i-th  radar  and  the  j-th  target.  (All  angles  are  defined  in  a 
horizontal  plane.) 

The  coordinates  describing  a given  configuration  will  be 
defined  as  components  of  a vector  X as  follows: 
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(the  angle  between  the  directions 
to  the  radars,  equal  to  the 
difference  of  the  radar  azimuths 
3^  with  respect  to  some  arbi- 
trary reference) 

x6  = ^1,1  = ^1  ~ ®i 

(©i  is  the  azimuth  of  target  i 
with  respect  to  the  arbitrary 
reference) 


= 02 


(Hi 


(the  altitude  of  OWN) 


is  the  altitude  of  target 


i) 


The  basic  BCAS  measurements  in  modified  form  are  the  com- 
ponents of  a vector  Y.  The  modifications  assure  that  the  mea- 
surement noise  in  each  term  is  independent  of  that  in  the  others 
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I 


and  has  variance  of  one.  It  may  be  observed  that  the  differen- 
tial azimuth  measurements  of  the  two  target  aircraft  for  the 
same  radar  have  correlated  errors.  Hence,  the  sums  and  the 
differences  of  the  differential  azimuths  are  used,  since  the 
errors  in  these  will  be  shown  to  be  uncorrelated.  The  proof  is 
the  following: 

Let 

5. 4.  Ill 


i.J 


0,  , - 0.  . 
i.J  i.o 


where  0,  , and  0,  are  the  bearings  of  the  target  and  own  air- 
l > J i . o 


craft  from  radar  i,  measured  with  respect  to  some  arbitrary 


reference.  a,  , is  in  fact  measured  as  the  difference  of  two 
1 > J 


such  angles.  The  error  in  measurement  for  each  is  the  error  in 
defining  the  radar  interrogator  beam  center. 

Hence,  the  measured  DAZ  is  5.4.15 


“ i.J  (0i,J  + ni,j)  " (0i,o  + ni,o) 


where  n.  . and  n.  are  zero-mean,  independent  noise  terms  as- 
i.J  i.o 


2 12  2 
sumed  to  have  equal  variances,  a = 2°a,  where  is  the  vari- 


ance of  the  DAZ  measurement  error,  estimated  from  flight  test 

data.  Then  the  noise  component  of  a.  . is  n = n.  , - n. 

i.J  “ij  i.J  i.o 


The  covariance  of  n 


and  n 


*i,l  ai,2 

E({nai,inai,i:  } )=  E({(ni,l-ni,o)(ni,2-ni,o)>  } 


is 


= E {ni,l  ni,2  “ ni,o  ni , 2 ~ ni,o  nil  + ni,o>  5.4.16 


= o^O 

since  the  noise  components  n^  , , n.  5 and  n.  are  mutually 

1)1  1 ) U 1 J O 


J 
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independent.  Hence  it  is  seen  that  the  measurement  noise  of 

a.  . is  correlated  with  that  of  a.  _. 

i , 1 1,2 


Consider  now  the  variables 


S = ai,l  + °i,2 


5-4.17 


^ ai,l_0li,2 


5.4.18 


<01,1  - 'll>0>t<‘Si,2  - «l.o) 
0i.l  * ei,2  - ®i,o  5.4.19 


and  the  noise  component  of  s is 


(3-4.2-10) 


which  has  the  variance 


ns  = ni,l  + ni ,2  " 2ni,o  5.4.20 


°s  = E 


IV] 


5.4.21 


Similarly, 


E jTn^  + ni>2  - 2ni}Q)2J 

E ti!l  + «i,2  + 4ni!o 
+2ni,lni12-J,ni,lnifo-i*ni,2ni,o  ] 


6 oP  = 3o  2 

a 


(0i,l^i,o)  * <0i,2  " 0i,o> 


= 0i,l  - 0i,2 

and  the  noise  component  of  d is  nd  with  a variance 

o2  = 2o2  = o2 
d a 


5.4.22 


5.4.23 
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■1,1 


5.4.25 


y6 


1,2 


l2,1 


y8 


2,2 


■55“  H» 


'10 


'11 


where  aa>  oT  and  oH  are  the  standard  deviations  of  the  DAZ,  TOA 
and  altitude  measurement  errors,  respectively. 

There  is,  then,  a set  of  non-linear  equations 


f±(X) 


5.4.26 


which  in  the  vicinity  of x =XQ  may  be  approximated  by  a set  of 
linear  equations. 


Ill 


5.  *<.27 


11  3f  (X) 

= W + * Tx  AXn 

n=l  n 


which  can  be  written  in  the  form 


11  3f,(X) 

A^i  = Z -37 Axn 

n=l  n 


5.4.28 


where  Ay  and  Ax  are  the  deviations  from  the  values  at  the  assumed 
configuration  X . The  set  of  equations  5.4.28  corresponds  in  form 
to  the  set  5.4.1  for  which  the  solution  was  derived  in  the  previous 
section.  The  partial  derivatives  evaluated  at  a.  point  near  the 
assumed  solution  play  the  roles  of  the  B's.  The  non-zero  partial 
derivatives  of  TOA,  DAZ,  and  altitude  with  respect  to  the  coordin- 
ates are  listed  in  Table  5.4-1. 

5.5  EFFECTIVENESS  OF  THE  INITIALIZATION  ALGORITHM 

It  may  be  noted  that  the  overall  process  of  finding  the  solu- 
tion consists  of  (1)  a coarse  initialization,  (2)  an  iterative 
improvement  on  the  initialization  and  (3)  an  iterative  squared 
error  minimization.  This  is  a complex  process,  and  one  is 
naturally  led  to  consider  whether  it  can  not  be  simplified. 

The  function  of  the  steps  is  the  following:  the  coarse 
initialization  (1)  is  necessary  to  establish  the  approximate 
region  in  the  space  of  configuration  coordinates  for  the  more 
refined  further  steps  to  be  applicable. 

The  iterative  squared  error  minimization  (3)  will  converge 
to  the  same  solution  independent  of  the  precise  initialization  as 
long  as  the  initial  point  is  within  a sufficiently  small  region 
about  the  solution  point  (a  local  minimum  of  the  summed  residuals). 
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TABLE  5.^-1.  PARTIAL  DERIVATIVES  OF  BCAS  MEASUREMENTS  WITH 
RESPECT  TO  COORDINATES  ABOUT  OWN 


TABLE  5-4-1.  PARTIAL  DERIVATIVES  OF  BCAS  MEASUREMENTS  WITH 
RESPECT  TO  COORDINATES  ABOUT  OWN  (Continued) 


aa1.1  - 

'"I  ^ 

COS*i.1  + r1 

30j 

4 - 2di 

rj  cos*ij  + 

4 

302,1  _ 

d2  r, 

2 

cos<Ji2^  - r^ 

3X  ~ 

dj  - 2d2 

r j cos*2j  + 

4 

Altitude  Measurement  Derivatives: 


(j  = 0,1,2) 


All  other  partial  derivatives  of  the  measurements  with  respect  to 
the  coordinates  are  identically  zero. 


If  the  initial  approximation  is  too  far  from  the  proper  solution, 
the  iteration  may  converge  to  an  entirely  different  solution, 
corresponding  to  some  other  local  minimum  of  the  summed  residuals. 
The  appropriate  measure  of  effectiveness  of  an  initialization 
scheme  in  this  context  then  is  not  the  accuracy  with  which  it 
itself  determines  the  configuration  of  the  aircraft  and  radars, 
but  rather  the  probability  that  the  least  square  minimization  which 
follows  will  converge  to  the  proper  solution.  That  solution  is  a 
maximum  likelihood  estimate  of  the  configuration,  given  the  measure- 
ments, and  is  independent  of  the  algorithm  used  to  find  it. 

The  coarse  initialization  step  (1)  was  not  expected  to  lead 
to  particularly  accurate  approximations  of  the  solution.  It  was 
for  this  reason  that  the  iterative  procedure  for  improving  the 
initialization  (2)  was  developed.  This  procedure,  too,  will  either 
converge  to  its  best  fit  to  the  data  if  it  is  started  close 
enough  to  it,  or  to  some  other  solution  otherwise.  The  regions 
of  convergence  of  the  two  iterative  procedures  are,  however, 
different.  Again,  if  this  procedure  to  improve  the  coarse  initial- 
ization is  used  in  conjunction  with  the  final  squared  error 
minimization  procedure,  the  measure  of  its  effectiveness  is  the 
probability  that  the  overall  process  ultimately  converges  to  the 
proper  solution. 

When  the  three  processes  are  applied  successively,  a number 
of  things  may  happen: 

a)  The  first  coarse  initialization  may  be  so  bad  that  no 
convergence  to  the  proper  solution  can  take  place  - i.e.,  the 
iterative  processes  will  converge  to  some  other  solution. 
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b)  The  first  coarse  initialization  may  not  give  a good 
starting  point  for  the  final  least-square  minimization  process, 
but  the  subsequent  iterative  improvement  to  the  initialization 
(2)  may  result  in  such  a "good"  starting  point. 

c)  The  iterative  squared  error  minimization  process  may 
converge  to  the  proper  solution  starting  from  either  the  approxi- 
mate solution  produced  by  process  (1)  or  the  presumably  improved 
approximation  of  process  (2).  It  may  take  fewer  iterations  to 
converge  from  one  than  from  the  other. 

d)  The  iterative  squared  error  minimization  may  converge  to 
the  proper  solution  if  started  at  the  approximation  produced  by 
the  first  coarse  iteration  (1),  but  not  if  started  at  the  approxi- 
mation produced  by  the  iterative  process  (2). 

Since  both  the  first,  coarse  initialization  (1)  and  the 
iterative  improvement  on  it  (2)  are  heuristic  algorithms,  it  was 
considered  necessary  to  test  them  experimentally  in  a range  of  con- 
figurations. The  performance  of  both  was  compared  to  what  could 
be  achieved  if  the  true  configuration  were  used  as  the  starting 
point  for  the  next  iterative  stage.  (The  true  configuration 
should  not  itself  be  the  solution  produced  by  any  of  the  algorithms, 
since  it  is  not  in  general  consistent  with  the  noise-corrupted 
measurements  from  which  the  algorithms  derive  their  solution.) 

The  following  series  of  computations  then  were  performed  for 
a range  of  configurations  and  the  results  were  compared: 

A.  (1)  Coarse  Initialization 

(2)  Improved  Initialization 

(3)  Residual  Minimization 
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B:  (1)  Coarse  Initialization 

(3)  Residual  Minimization 

C:  (0)  Initialization  with  True  Configuration 

(3)  Residual  Minimization 

D:  (0)  Initialization  with  True  Configuration 

(2)  Improved  Initialization 

(3)  Residual  Minimization 

It  is  also  of  interest  to  see  how  well  the  initialization 
algorithms  (1)  and  (2)  approximate  the  true  configuration,  since 
if  their  performance  were  reliably  good,  the  final  step  (3)  might 
be  an  unnecessary  refinement. 

The  configurations  for  which  these  tests  were  performed  are 
shown  in  Figures  5.5-1  through  5.5-3-  In  each  case,  the  BCAS  air- 
craft and  one  intruder  aircraft  were  assumed  to  be  at  fixed  posi- 
tions relative  to  a configuration  of  two  radars.  Results  were 
then  evaluated  with  a second  intruder  aircraft  placed  in  a sequence 
of  positions  on  a circle  about  the  BCAS  aircraft.  At  each  posi- 
tion, the  calculation  was  performed  ten  times,  with  a different 
set  of  simulated  random  measurement  errors  each  time.  A represent- 
ative set  of  outputs  for  one  configuration  is  shown  in  Figures 
5.5-^a  to  5-5-^c. 

Qualitatively,  the  results  were  the  following:  Im  most  con- 
figurations, the  final  squared  error  minimization  process  converges 
to  the  proper  solution,  independent  of  what  process  is  used  to 
arrive  at  the  starting  approximation. 

There  are  some  configurations  In  which  neither  sequence  A nor 
the  sequence  B for  obtaining  the  starting  approximations  allowed 
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138.59 


Hq  = 20,000' 


FIGURE  5-5-2.  SIMULATED  PATTERN  2. 
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FIGURE  r;.5-Jlb.  SAMPLE  OUTPUT  OF  SIMULATION  RUN  (CONTINUED). 


FIGURE  5.5-4c.  SAMPLE  OUTPUT  OF  SIMULATION  RUN  (CONTINUED). 


the  final  process  to  converge  to  a practically  useful  arrroxi mation 
of  the  actual  configuration. 

Overall,  the  final  process  is  most  likely  to  converse  when  it 
begins  at  the  approximate  solution  ierived  hy  the  initial  "oarse 
approximation  followed  hy  the  iterative  improvement  (sequence  A). 
However,  the  si*  nation  dees  ,-ccur  that  th---  final  squared-error 
nini r**.i 2at i on  mocedure  conv^r’^r  to  the  merer  solution  when 
starting  from,  the  initial  "coarse"  approximation,  but  not  from  the 
configuration  obtained  by  *-he  iterative  "improvement"  (sequence  E). 

Solutions  were  attempted  a*:  216  different  configurations  of 
radars  and  targets.  Ten  sets  of  noisy  measurements  were  simulated 
at  each,  position.  In  106  of  these  con.f  i gurat  ions , the  final 
squared  error  minimisation  process  converged  to  the  proper  solution 
every  time,  regardless  of  whether  the  initial  coarse  approximation 
cr  the  iteratively  improved  approximation  were  used  as  starting 
points.  In  159  conf imurat ions , the  rrorer  solution  was  obtained 
every  time  if  the  better  of  the  two  rcssitle  initializations  were 
selected . 

If  the  success  rate  is  calculated  over  the  total  number  of 
trials  - 216  positions  with  10  sets  of  simulated  noisy  measurements 
at  each  - the  following  percentages  are  obtained: 

Successful  Convergence  with  either  Initialization 
■ A and  E)  69.9? 

_ Successful  Convergence  onlv  v 1th.  Fins!  7r.lt  isll  zsit  icr 
(3,  not  A)  r.P5' 

Successful  Convergence  Only  with  Iteratively  Improved 

Initialization  (A,  not  E)  7.6? 
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No  Convergence  with  Either  Initialization  (Neither 


A nor  B)  17*7? 

The  squared  error  minimization  process  (3)  converged  when  the 
actual  configuration  was  used  as  the  initialization  more  than  with 
any  other  initialization  (though  this  does  not  imply  that  the 
solution  found  was  always  a good  approximation  to  the  actual  con- 
figuration). There  were  also  configurations  in  which  the  sequences 
involving  step  (2),  the  iterative  improvement  of  the  initialization 
(i.e.,  sequences  A and  D)  diverged,  whereas  all  other  sequences 
converged.  These  observations  indicate  that  the  heuristically 
derived  initialization  algorithms  could  possibly  be  Improved.  On 
the  whole,  however,  the  initialization  algorithms  tend  to  fail  in 
what  are  inherently  bad  ranges  for  the  system  of  equations  - that 
is,  configurations  in  which  the  situation  is  inherently  highly 
sensitive  to  measurement  errors.  Hence,  improvement  in  the  initial- 
ization procedures  would  not  necessarily  lead  to  significantly 
improved  system  performance. 

The  approximate  configurations  computed  by  the  initialization 
schemes  on  the  average  are  considerably  farther  from  the  true 
configurations  than  the  configurations  computed  by  the  final  step 
which  minimizes  the  residual.  Hence,  the  initialization  algorithms 
by  themselves  do  not  give  results  indicative  of  how  well  the  con- 
figuration can  be  determined  from  the  measurements.  The  values 
that  they,  by  themselves,  determine  should  not  be  used  for  target 
tracking. 
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5.6  OVERALL  ASSESSMENT  OF  THE  SYSTEM  WITHOUT  AZIMUTH  REFERENCE 

SIGNALS 

The  accuracy  of  the  target  range  and  bearing  calculations  that 
can  be  obtained  from  the  BCAS  measurements  based  on  radars  without 
azimuth  reference  signals  can  be  judged  from  the  results  of  the 
simulations  shown  in  Section  6. 

As  in  the  case  of  radars  with  azimuth  reference  signals, 
there  are  ranges  of  configurations  in  which  good  solutions  can  not 
be  obtained.  These  configurations  include  the  following: 

(1)  when  the  two  radars  are  colinear  when  viewed  from  the 
BCAS 

(2)  when  either  of  the  target  aircraft  is  in  between  the 
BCAS  and  one  of  the  radars 

(3)  when  the  BCAS  is  directly  between  either  target  aircraft 
and  one  of  the  radars 

(4)  when  the  BCAS  aircraft  and  both  target  aircraft  are  in  a 
line . 

The  extent  of  each  band  range  is  a function  of  the  total  con- 
figuration. 
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6.  SIMULATION 

6.1  SIMULATION  TECHNIQUE 

The  various  algorithms  for  determining  radar  and  target  range 
and  bearing  from  own  on  the  basis  of  the  BCAS  measurements  were 
tested  in  a series  of  digital  computer  simulations. 

In  principle,  all  the  simulations  contained  the  following 
steps.  A configuration  of  radars  and  own  and  target  aircraft  was 
assumed.  The  TOA's,  DAZ ' s , and  OAZ’s  that  the  BCAS  would  measure 
in  such  a configuration  in  the  absence  of  noise  and  other  sources 
of  error  (such  as  quantization  effects)  were  computed.  Random 
numbers  were  generated  to  represent  the  various  error  sources  and 
were  added  to  the  "perfect  measurements".  The  resulting  "noise- 
corrupted  measurements"  were  then  used  as  inputs  to  the  computa- 
tional algorithms,  which  used  them  to  derive  the  positions  of  the 
radars  and  targets  relative  to  OWN.  The  results  of  these  computa- 
tions were  compared  with  the  originally  assumed  configuration  to 
find  the  errors. 

The  random  error  terms  added  to  the  measurements  were  drawn 
from  normal  distributions  with  variances  chosen  to  match  the 
variances  of  the  corresponding  measurements,  as  determined  during 
the  test  flights.  Where  angular  measurements  have  correlated 
errors,  the  random  error  terms  in  the  simulations  reflected  this. 
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The  errors  introduced  were  the  following: 

TOA:  Normally  distributed  error,  mean  zero,  standard 

deviation  0.15  microseconds,  independent  for  each  measurement. 

Altitudes:  Normally  distributed  error,  mean  zero, 

standard  deviation  30  ft.,  independent  for  each  measurement. 
(This  is  probably  not  a good  way  to  model  altitude  errors; 
however,  the  computed  results  were  quite  insensitive  to  alti- 
tude errors,  so  that  the  faults  of  the  model  are  not  important.) 

Azimuth  measurements:  The  measurement  error  is  due  to 
the  incorrect  determination  of  the  beam  center  position  from 
the  bursts  of  pulses  received  when  the  beam  is  pointed  at  OWN, 
at  OTHER,  or  in  an  azimuth  reference  direction.  The  errors 
were  modeled  as  zero-mean,  normally  distributed,  random  compo- 
nents with  standard  deviation  of  0.2°  in  each  measurement 
of  direction.  This  leads  to  DAZ  and  OAZ  errors  with  standard 
deviations  of  0.3°  and  the  proper  mutual  correlation. 

Although  the  simulations  were  set  up  in  such  a way  that 
one  or  more  of  the  aircraft  were  placed  at  consecutive  posi- 
tions of  what  might  be  flight  trajectories,  the  calculations 
at  each  position  were  carried  out  completely  independent  of 
each  other.  The  purpose  of  the  simulations  was  to  determine 
how  well  the  BCAS  measurements  at  a single  point  could  deter- 
mine the  aircraft  configuration.  An  actual  operational  system 
should  do  better  since  it  would  smooth  the  data,  i.e.,  combine 
the  measurements  at  a given  instant  with  previous  measurements 
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in  such  a way  as  to  give  a better  estimate  of  the  true  posi- 
tion. The  simulation  in  each  case  used  the  minimum  number  of 
measurements  to  determine  the  configuration.  Again,  an  opera- 
tional system  might  frequently  do  better,  using  more  than  the 
minimum  number  of  radars  to  track  a target  aircraft,  or  basing 
calculations  of  radar  position  on  more  than  the  minimum  number 
of  tracked  targets. 

Since  the  simulations  were  strictly  for  static  position 
calculations  rather  than  for  dynamic  computation  of  aircraft 
tracks,  one  complicating  factor  with  regard  to  the  measurements 
was  neglected.  It  was  assumed  that  the  measurements,  subject 
to  the  errors  discussed,  were  available  at  the  nominal  times. 

In  the  actual  BCAS  environment,  the  measurements  are  made  at 
times  determined  by  the  asynchronous  rotation  of  the  locked 
radars,  i.e.,  roughly  every  four  seconds  for  the  short-range 
SSR's.  Thus,  when  the  measurements  from  one  radar  are  re- 
ceived, those  from  the  other  may  be  up  to  about  four  seconds 
old.  An  operational  system  would  have  to  include  some  mecha- 
nism for  taking  the  differences  in  measurement  times  into 
account.  Experience  with  test  data  (see  Section  7)  indicates 
that  the  effect  of  the  asynchronous  measurement  times  is 
minor.  Therefore,  the  results  of  the  static  simulation  should 
not  be  significantly  affected  by  omitting  this  factor  from 
consideration. 
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6.2  THE  SITUATIONS  SIMULATED 

The  purpose  of  the  simulations  was  to  test  quantitatively 
the  accuracy  of  the  BCAS  range-bearing  calculating  algorithms 
for  the  cases  with  and  without  SSR  azimuth  reference  signals. 

A succession  of  configurations  with  two  radars,  a BCAS 
aircraft,  and  two  target  aircraft  were  set  up  and  the  "noise- 
free  measurements"  of  TOA,  DAZ  and  OAZ  were  calculated.  To 
this  set  of  "measurements"  appropriate  "noise"  terms  were 
added  to  create  a set  of  "noisy  measurements".  Each  algorithm 
then  used  the  appropriate  subset  of  these  "noisy  measurements" 
to  compute  the  range  and  bearing  from  OWN  of  the  target  air- 
craft and  the  radars.  The  errors  of  the  resulting  solutions 
were  computed  by  comparing  the  solutions  obtained  with  the 
configuration  initially  assumed. 

The  calculations  using  azimuth  reference  pulses  from  two 
radars  were  carried  out  for  each  target  aircraft 
separately.  The  calculations  using  the  "no  azimuth  reference" 
solution  were  carried  out  for  the  configuration  as  a whole. 

Two  different  situations  were  assumed  for  this  case.  These 
were  (1)  partial  equipage  of  radars  with  azimuth  reference 
signals  and  (2)  no  azimuth  reference  signals  from  any  ground 
radars. 

For  the  first  situation,  it . is  assumed  that  in  each  radar 
coverage  region,  one  SSR  is  equipped  with  azimuth  reference 
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signals.  Then,  once  the  BCAS  has  computed  the  relatiye  direc- 
tions of  the  radars  and  target  aircraft,  it  can  use  the  di- 
rectly measured  direction  to  this  one  radar  to  properly  orient 
the  whole  configuration. 

For  the  second  situation,  it  is  assumed  that  the  BCAS 
computes  the  shape  of  the  configuration,  but  has  no  direct  way 
of  establishing  its  orientation.  This  can  be  done  only  in  a 
dynamic  situation  when  the  BCAS  is  flying  in  a known  direction 
and  computing  the  shape  of  the  configuration  at  successive 
points  along  its  path.  The  distance  between  the  two  SSR's  and 
the  direction  from  one  SSR  to  the  other  do  not  change  with 
time.  In  concept,  the  BCAS  can  then  determine  its  own  position 
relative  to  this  fixed  line  at  a number  of  successive  points 
in  time.  The  angle  between  the  fixed  line  between  the  radars 
and  the  line  determined  by  the  successive  positions  of  the  BCAS 
aircraft  (which  lies  along  the  known  flight  direction)  can  then 
be  determined  and  used  to  orient  the  whole  configuration  of 
radars  and  aircraft. 

For  the  first  of  these  two  situations,  the  measures  of 
accuracy  for  the  computed  results  are  the  RMS  errors  in  range 
and  bearing  to  the  radars  and  the  other  aircraft.  For  the 
second  situation,  the  range  errors  are  inherently  the  same  as 
for  the  first.  The  bearing  errors  are  computed  relative  to  the 
line  connecting  the  radars,  which  is  assumed  known. 
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6.3  SIMULATION  RESULTS 


The  patterns  of  simulated  configurations  of  radars  and  air- 
craft are  shown  in  Figures  6.3-1  to  6.3-^.  The  aircraft  were  placed 
at  successive  positions  of  their  "tracks"  as  though  they  were 
moving  in  straight  lines  at  a constant  speed.  In  all  but  the 
final  pattern,  there  was  a "near  miss"  between  the  BCAS  and 
Target  1. 

Figures  6.3-5  to  6.3-38  show  the  errors  in  the  computed  range 
and  bearing  values,  that  is,  the  differences  between  the  com- 
puted and  true  values.  All  points  that  would  have  been  outside 
the  limits  of  the  plots  (in  either  direction)  are  shown  below 
the  horizontal  axis.  These  may  be  regarded  as  "wild"  points 
in  the  results. 

All  the  values  of  a given  variable  which  were  computed  by 
different  techniques  are  shown  on  the  same  plot,  indicated  by 
different  symbols.  It  is  to  be  noted  that  when  radar  distances 
were  computed  using  the  azimuth  reference  signals  from  both 
radars,  a separate  calculation  was  carried  out  for  each  target. 

The  results  are  shown  separately.  The  radar  bearings  for  this 
case  are  then  merely  the  simulated  measurement  errors. 

The  bearing  errors  shown  for  the  computations  based  on  no 
azimuth  reference  signals  are  calculated  relative  to  the  line 
connecting  the  radars. 

When  one  radar  only  is  assumed  to  have  azimuth  reference 
signals,  it  is  always  radar  1. 
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Pattern  1 

For  Pattern  1 (Figure  6.3-1  and  Figures  6.3-5  to  6.3-12),  the 
initial  configurations  allow  accurate  computations  of  most 
variables  by  all  the  computational  schemes.  In  some  cases,  the 
no-azimuth-reference  algorithm  fails  to  converge  because  of  bad 
initialization.  This  occurs  because  the  BCAS  is  almost  di- 
rectly between  radar  2 and  target  1.  Where  the  algorithm  over- 
comes this  difficulty,  the  computed  results  are  very  accurate. 

The  final  configurations  are  generally  unfavorable  to 
accurate  computation,  largely  because  target  1 comes  between 
the  BCAS  and  both  target  2 and  radar  2. 

The  computations  of  distance  to  radar  2 are  bad  throughout 
because  the  differential  azimuths  of  both  targets  with  respect 
to  radar  2 are  consistently  small. 

The  bearing  errors  are  seen  to  be  large  for  the  no-azimuth- 
reference  system  and  quite  acceptable  with  two  radars  with  azi- 
muth references.  The  results  for  the  system  with  one  azimuth 
reference  signal  are  intermediate. 

Pattern  2 

For  Pattern  2,  good  solutions  are  found  by  all  systems  in 
all  but  three  limited  conditions. 

The  initial  configuration  is  such  that  radar  1 is  very 
close  to  the  BCAS  and  between  the  BCAS  and  target  1.  BCAS 
operation  in  such  a configuration  is  not  anticipated.  All 
systems  fail  to  compute  the  configuration. 

133 


J 


IGURE  6.3 
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When  the  BCAS  aircraft  and  target  1 are  in  immediate 
proximity  of  each  other,  the  solutions  fail.  Finally,  at  both 


ends  of  the  pattern  there  appears  to  be  a false  solution  for 
the  position  of  target  2 computed  using  the  technique  utilizing 
only  that  target  and  azimuth  references  on  both  radars. 

Where  good  solutions  are  obtained,  the  accuracy  of  the 
computed  distance  to  radar  1 (the  nearby  radar)  is  very  good 
for  both  systems  of  computation.  The  accuracy  of  the  computed 
range  to  radar  2 (the  far  radar)  is  much  worse.  It  is  better 
for  the  no-azimuth-reference  system  with  two  targets  than  for 
the  system  with  azimuth  references  applied  to  each  target 
separately . 

The  errors  in  the  computed  target  distance  are  comparable 
for  the  two  systems  and  are  mostly  less  than  200  feet. 

Target -bearing  errors  are  comparable  for  all  the  systems, 
with  an  RMS  value  of  less  than  one  degree  (except  when  the 
target  is  quite  close).  Radar-bearing  errors  (for  the  systems 
that  do  not  measure  OAZ  directly)  are  larger.  In  this  instance, 
the  error  in  bearing  to  radar  2 is  greater  if  the  OAZ  to  radar 
1 is  measured,  (one  azimuth  reference  case)  than  if  no  azimuth 
reference  is  measured,  but  the  radars  are  assumed  to  lie  in  a 
known  direction,  one  from  the  other.  This  atypical  result  is 
due  to  the  very  high  accuracy  with  which  the  distance  to  radar 
1 is  computed,  and  the  insensitivity  of  the  angle  computations 
to  the  distance  to  radar  2. 
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For  the  configurations  that  constitute  pattern  3,  the  algor- 
ithm depending  on  azimuth  reference  signals  from  both  radars  con- 
verges to  a satisfactory  solution  almost  everywhere.  The  exceptions 
are  two  individual  points  near  the  beginning,  when  the  calculations 
involving  target  1 fail  to  converge.  It  must  be  noted  that  in 
these  instances  BCAS,  the  target,  and  radar  1 are  at  the  vertices 
of  an  approximately  equilateral  triangle  - a configuration  that 
violates  the  assumptions  about  the  nature  of  the  solutions  that  are 
to  be  sought. 

It  may  also  be  noted  that  at  about  the  point  of  nearest 
approach  (at  the  time  of  the  near  miss),  the  computed  target  bear- 
ing and  radar  distance  values  have  large  errors,  but  the  target 
distance  does  not. 

In  contrast,  it  may  be  noted  that  the  no-azimuth-reference 
solutions  for  the  configurations  near  the  beginning  of  the  pattern 
show  large  errors  in  the  computed  distance  to  target  1,  while  the 
other  computed  quantities  do  not  have  exceptionally  large  errors. 
This  example  illustrates  the  rather  subtle  nature  of  the  "bad 
ranges"  that  may  be  encountered  in  solving  for  the  shape  of  the 
configuration  in  the  absence  of  azimuth  reference  signals.  All 
the  solutions  obtained  are  correct  in  the  sense  that  the  computed 
residuals  show  that  they  fit  the  measurements  quite  closely.  The 
great  error  sensitivity  seems  to  arise  from  the  circumstance  that, 
when  viewed  from  the  BCAS,  both  targets  (in  direction)  lie  between 
the  radars,  33°  and  H6°  from  the  radars.  When  these  angles  become 
slightly  larger,  the  error  sensitivity  of  the  solution  disappears. 
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FIGURE  6.3-3.  SIMULATED  PATTERN  3. 
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The  range  errors  to  both  targets  become  approximately  200  feet 
RMS. 


Pattern  4 

' 

Pattern  4 illustrates  the  case  of  an  alternate  solution  being 
found  by  the  technique  relying  on  both  radars  being  equipped  with 
azimuth  reference  signals.  The  false  solution  is  consistently 
found  for  target  2 from  the  34th  sampling  point  on.  The  smoothness 
and  continuity  of  the  false  track  is  evident  in  the  error  plot  for 
the  computed  distance  to  radar  1 (Figure  6.3-29),  and  the  distance 
and  bearing  to  the  target  on  the  plots  for  which  the  scale  has 
been  compressed  (Figures  6.3-33  and  6.3-38). 

The  other  effects  seen  for  this  pattern  that  merit  comment 

are : 

The  differential  azimuth  measurements  relative  to  radar  2 are 
small  for  both  targets.  Correspondingly,  the  errors  in  the  com- 
puted  distance  to  radar  2 are  large.  This  has  the  further  effect 
of  making  the  error  in  the  computed  bearing  of  this  radar  (for  the 
no-azimuth-reference  and  one-azimuth-reference  cases)  large. 

For  the  later  configurations  in  the  pattern,  both  targets  are 
approximately  between  the  BCAS  and  radar  1.  In  these  configura- 
tions, the  no-azimuth-reference  solution  has  very  large  errors  in 
computed  range.  The  bearing  errors  appear  to  have  a somewhat 
narrower  bad  range. 

In  the  same  set  of  configurations,  the  2-azimuth-reference 
solution  for  target  1 is  also  bad,  again  because  target  1 is 
between  the  BCAS  and  the  radar. 
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FIGURE  6.3-14 
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FIGURE  6.3-22. 


FIGURE  6.3-23. 
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7.  FLIGHT  TEST  DATA 


Representative  riots  of  ranre  and  bearing  to  a stationary 

a 

target  compute  i from  actual  measurements  made  In  flight  are  shovm 
in  Figures  7.1-1  and  7.1-2.  These  represent  calculated  range  and 

1 S 

bearing  to  the  Mizpah  fire  tower  from  the  PCAS  test  aircraft 
flying  over  it.  The  calculations  are  based  on  measurements 

utilizing  two  radars  v/ith  azimuth  reference  pulses.  The  circles  jj 

represent  active  interrogator  measurements  of  range,  the  x's  the 

range  and  bearing  computed  by  BCAS,  the  short  horizontal  bars  the 

range  and  bearing  computed  from  ARTS  III  observations,  and  the 

vertical  bars  through  them  the  approximate  90*  confidence  limits 

for  the  values  computed  from  the  ARTS  measurements. 


FIGURE  7.1-1.  COMPUTED  RANGE  TO  TARGET. 


TEST  DATE:  OCT.  15,  1975 
RUN  #1  PATTERN  #10  MODE  133 


SUMMARY 


Analyses  have  been  performed  for  the  static  case  and  algor- 
ithms have  been  developed  for  computing  range  and  bearing  to 
potential  threat  aircraft  and  to  ground  radars,  using  passive-mode 
BCAS  measurements  as  well  as  a combination  of  active  and  passive 
mode  measurements  for  the  single  ground  site  geometry.  Only  the 
static  solutions  are  presented.  These  are  solutions  based  on  the 
differential  time  of  arrival  (TOA),  differential  azimuth  (DAZ) 
and,  where  appropriate,  own  azimuth  (OAZ)  measurements  taken  at 
one  instant  for  the  configuration  as  it  exists  at  that  time. 

The  algorithms  compute  range  and  bearing  to  intruder  aircraft 
and  radars  on  the  basis  of  only  the  current  measurements,  making 
use  of  no  a priori  knowledge  of  either  the  positions  of  the  air- 
craft or  radars  or  of  any  previous  measurements.  For  this  reason, 
the  accuracies  of  the  computed  positions  are  likely  to  be  worse 
than  those  that  would  be  obtained  by  dynamic  tracking  algorithms 
which  would  smooth  out' the  effects  of  measurement  errors  over 
time . 

The  solutions  obtained  are  the  solutions  that  best  fit  the 
measurement  data.  Thus,  their  accuracy  is  the  intrinsic  accuracy 
to  within  which  the  relative  positions  of  intruder  aircraft  and 
locked  radars  can  be  determined  from  the  measurements  made  with  the 
given  accuracy  at  one  point  in  time.  Dynamic  tracking  algorithms, 
when  they  are  developed,  may  be  expected  to  have  better  overall  per- 
formance because  they  will  have  available  sequences  of  positions 
over  a period  of  time  and  will  be  able  to  smooth  out  measurement 
errors  by  effectively  averaging  over  time.  The  statically  computed 
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accuracies  should  be  suggestive  of  the  accuracies  that  the  tracking 
algorithms  should  be  able  to  achieve. 

Algorithms  to  determine  target  range  and  bearing  from  the 
BCAS  measurements  were  developed  and  tested  in  simulations.  In 
addition,  an  approach  to  the  problem  using  target  locus  curves  has 
been  developed  both  as  an  aid  to  assessing  the  nature  of  the 
solutions  qualitatively  and  as  a potential  basis  for  improved 
solution  algorithms. 

Three  different  modes  of  purely  passive  operation  have  been 
simulated.  These  assume  all  radars  equipped  with  azimuth  refer- 
ence signals,  none  so  equipped,  and  only  one  radar  so  equipped 
available  at  a given  time. 

It  appears  that  when  no  radars  are  equipped  with  azimuth 
references,  the  target  bearing  cannot  be  derived  sufficiently 
accurately  to  give  good  target  tracks.  This  conclusion  should  be 
verified  by  dynamic  simulation. 

When  all  radars  are  equipped  with  azimuth  reference  signals, 
the  positions  of  single  targets  can  be  determined.  The  range  of 
configurations  in  which  the  solution  is  excessively  error- 
sensitive  is  smaller  than  for  the  other  cases,  but  under  some 
circumstances  the  measurements  lead  to  ambiguities  in  that  several 
distinct  configurations  can  give  rise  to  the  same  measurements. 

When  only  one  radar  has  azimuth  reference  signals,  two  target 
aircraft  must  be  observed  to  make  calculations  of  position  based 
on  passive  measurements  possible.  The  range  of  configurations  in 
which  the  solution  is  excessively  sensitive  to  measurement  error 
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Is  larger  than  when  all  radars  have  azimuth  reference  signals; 
multiple  solutions  do  not  occur. 


In  the  good  configurations  of  radars  and  aircraft,  target 
aircraft  positions  can  be  determined  to  within  an  RMS  error  of 
less  than  300  feet,  assuming  measurement  accuracies  like  those 
obtained  by  the  experimental  BCAS  system.  Radar  positions  can  be 
determined  only  to  within  a few  miles  using  the  system  of  two 
radars  with  azimuth  reference  and  one  target;  radar  position  error 
is  under  one  mile  when  the  system  with  two  target  aircraft  is 
used. 

It  is  judged  that  either  system  - assuming  all  radars  equipped 
with  azimuth  reference  signals  or  only  one  radar  within  BCAS  range 
so  equipped  - Is  technically  feasible. 

The  results  for  the  individual  cases  are  summarized  below: 

1.  All  radars  are  assumed  equipped  with  azimuth  reference 
signals . 

In  this  situation,  the  range  and  bearing  from  BCAS  of  a 
single  Intruder  aircraft  can  be  determined  if  It  is  being  tracked 
by  BCAS  using  two  or  more  ground  radars.  The  range  to  each  radar 
can  also  be  calculated. 

Both  simulated  and  flight  test  data  verify  that  this  mode  of 
operation  is  possible  in  all  but  a set  of  unfavorable  configura- 
tions. The  unfavorable  configurations  are  the  following: 

a)  BCAS  on  an  extension  of  the  line  connecting  the  two  radars. 

b)  The  intruder  aircraft  between  BCAS  and  either  of  the 
radars . 

c)  The  BCAS  aircraft  between  one  of  the  radars  and  the 
intruder. 
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The  width  of  the  bad  ranges  depends  on  the  characteristics 
of  the  configuration  as  a whole.  In  the  worst  part  of  each, 
the  iterative  solution  algorithm  fails  to  converge  to  any  solu- 
tion. Near  the  edges  of  the  bad  region,  the  configuration  com- 
puted from  the  measurements  is  highly  sensitive  to  measurement 
errors.  This  error  sensitivity  is  intrinsic  In  that  the  con- 
figurations are  such  that  large  changes  in  the  relative  posi- 
tions of  the  aircraft  (and  radars)  cause  small  changes  in  the 
observed  measurements.  Then,  inversely,  small  changes  in  the 
measurements,  such  as  those  arising  from  measurement  noise, 
cause  large  changes  in  the  configuration  that  can  be  deduced 
from  the  measurements. 

Outside  the  bad  ranges,  the  relative  position  of  the  in- 
truder aircraft  can  be  calculated  with  an  RMS  error  in  position 
of  generally  less  than  300  feet,  depending  on  the  configuration. 

The  ranges  to  the  radars  can  also  be  calculated,  but  the 
values  obtained  are  quite  sensitive  to  errors  in  the  measurement 
of  differential  azimuth  and  may  have  errors  of  several  miles. 

It  is  important  to  notice  that,  under  some  circumstances, 
several  distinct  configurations  of  radars  and  aircraft  will  produce 
the  same  set  of  values  for  all  the  measurements  obtained  by  the 
BCAS . In  such  a case,  it  is  theoretically  impossible  to  deter- 
mine from  thfe  set  of  static  measurements  obtained  at  one  time 
which  of  the  possible  configurations  actually  gave  rise  to  the 
measurements.  The  ambiguity  can  be  resolved  by  making  other 
measurements,  e.g.,  an  active  measurement  of  target  range.  In 
addition,  although  the  possibility  of  the  false  solution  may 


182 


persist  for  a period  of  time,  so  that  a false  track  for  the  in- 
truder may  be  established  instead  of  the  true  one,  the  radar 
positions  computed  in  conjunction  with  the  false  track  will 
in  time  be  seen  to  be  inconsistent  in  that  the  radars  will 
appear  to  move  relative  to  each  other. 

2.  No  radars  are  assumed  equipped  with  azimuth  reference 
signals . 

In  this  situation,  BCAS  can  determine  the  shape  of  the 
configuration  of  radars  and  aircraft  if  there  are  two  radars 
and  two  or  more  target  aircraft.  The  measurements  contain  no 
absolute  azimuth  reference  signal.  Hence,  the  orientation  of 
the  configuration  cannot  be  determined  directly,  but  only  the 
bearing  of  each  radar  and  aircraft  relative  to  some  arbitrary 
reference  within  the  configuration . 

It  has  been  suggested  that  the  BCAS-equipped  aircraft  might 
be  able  to  compute  its  own  position  relative  to  the  radars  at 
a number  of  consecutive  times.  Then  it  could  relate  its  own 
flight  direction  to  the  fixed  direction  of  the  line  connecting 
the  two  radars  and  use  that  as  the  known  reference  direction  in 
determining  the  bearing  angles  toward  the  intruder  aircraft. 

It  was  found  in  the  course  of  the  simulations  that  the 
range  of  configurations  in  which  the  computed  results  are  in- 
trinsically highly  sensitive  to  measurement  error  is  more  exten- 
sive in  this  case  than  in  the  case  where  both  radars  are 
equipped  with  azimuth  reference  signals.  The  bad  configurations 
include  the  following: 
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a)  BCAS  on  extension  of  line  connecting  the  two  radars. 

b)  Either  Intruder  aircraft  between  BCAS  and  either  of 
the  radars. 

c)  Both  intruder  aircraft  in  the  same  direction  as  viewed 
from  BCAS. 

In  addition,  the  current  heuristic  algorithm  used  to  obtain 
an  Initial  approximation  tends  to  fail  when  an  intruder  air- 
craft is  in  the  direction  directly  opposite  that  of  a radar, 
when  viewed  from  BCAS,  or  when  BCAS  is  directly  between  the 
intruder  aircraft.  The  difficulty  in  these  regions  can  be 
eliminated,  since  it  arises  from  the  algorithm  used  and  not  from 
the  nature  of  the  dependence  between  the  configuration  and  the 
measurements. 

For  those  configurations  in  which  the  aircraft  and  radar 
positions  can  be  reliably  computed,  (i.e.,  in  the  good  ranges), 
the  errors  in  range  to  the  target  aircraft  are  comparable  to 
the  errors  obtained  using  radars  that  are  all  equipped  with 
azimuth  reference  signals.  The  computed  radar  distances  are 
considerably  more  accurate.  However,  the  bearing  angles  com- 
puted relative  to  the  line  connecting  the  radars  have  errors  on 
the  order  of  several  degrees.  This  suggests  that  the  proposed 
scheme  of  operation  with  no  azimuth  reference  signals  at  all 
may  have  difficulties.  A definitive  judgment  must  rest  on  an 
analysis  of  a tracking  scheme  in  the  dynamic  situation. 

3.  Some  radars,  but  not  all,  are  equipped  with  azimuth 


reference  signals. 


It  is  assumed  that  the  BCAS  at  any  one  time  would  be  in 


range  of  only  one  radar  with  azimuth  reference  signals.  Other 
radars  would  be  available  for  locking  and  for  tracking  targets, 
but  these  would  not  have  azimuth  reference  signals. 

In  such  a situation,  the  problem  of  computing  the  configu- 
ration of  radars  and  aircraft  at  any  one  instant  is  essentially 
the  same  as  in  the  case  of  no  azimuth  reference  signals.  Two 
radars  and  two  aircraft  are  required.  The  only  difference  is 
that,  once  the  shape  of  the  configuration  has  been  determined, 
it  can  be  properly  oriented  on  the  basis  of  the  azimuth  measure- 
ment . 

The  extent  of  the  bad  range  for  this  case  is  identically 
the  same  as  in  the  case  of  the  system  with  no  azimuth  refer- 
ences, as  is  the  error  sensitivity  of  the  computed  ranges  to  the 
intruder  aircraft  and  the  radars.  The  bearing  errors  to  the 
radars  and  the  aircraft  are  comparable  to  those  achieved  when 
all  radars  are  equipped  with  azimuth  reference  signals. 

. Limited  simulations  were  also  performed  for  the  case  of 
a single-site  geometry,  where  active  interrogation  is  used  to 
establish  target  range.  It  was  found  that  the  accuracy  to  be 
obtained  is  comparable  to  that  for  the  purely  passive  cases. 

The  worst  accuracies  occur  when  the  radar,  the  BCAS,  and  the 
target  are  approximately  colinear.  Two  distinct  solutions  can 
occur  if  the  target  is  at  a sufficiently  high  altitude. 
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ALGEBRAIC  CONDITIONS  FOR  PROPERTIES  OF  LOCUS  CURVE 

A . 1 GEOMETRIC  ARGUMENTS 

The  properties  of  the  target  locus  curves  described  in  Section 
2 largely  depend  on  the  relative  values  of  T,  H and  HQ . The  single 
most  important  determinant  of  the  nature  of  the  locus  curves,  in 
terms  of  the  geometric  construction  carried  out  previously,  is 
whether  the  vertical  projection  of  the  radar  position  is  inside  or 
outside  the  ellipse  in  the  target  altitude  plane. 

Figure  A-l  shows  a side  view  of  the  ellipsoid  of  revolution. 


I 


cut  along  the  LOP.  The  critical  question  is  whether  x is  positive 
(as  shoyn)  or  negative.  For  positive  x,  the  radar  projection  is 
inside  the  ellipse;  for  x negative,  it  is  outside. 

The  quantities  in  the  figure  are  defined  by 


FIGURE  A-l.  SIDE  VIEW  OF  ELLIPSOID  OF  REVOLUTION  USED  TO 
CONSTRUCT  TARGET  LOCUS  CURVES. 
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s1  = yx2  + h2 

A1.2 

R1  = ^(d+x)2  + (H-Hq)2 

A1.3 

By  definition 

T = S1  + R1  - D 

Al.l» 

Hence 

(T+D-S1)2  = (R1)2 

A1.5 

or 

(T+D)2 

- 2 (T+D ) ^x2  + H2  + x2  + H2 

= d2  + 2dx  + x2  + H2  - 2HH  + H2 

0 o 

A1.6 

Hence 

T2  + 2TD  + 2HH 

0 

- 2dx  = 2 (T+D)  y^x2  + H2 

A1.7 

and,  after  squaring 

(T2+2TD+2HH  )2  ■ 

O 

- b (T2+2TD+2HH  ) dx 
o 

+ lJd2  x2 

= 4 (T+D ) 2 (x2+H2) 

A1.8 

or 

^ 4 (T+D ) 2 - i)d2] 

x2  + b ( T2+2TD+2HHq ) dx 

+ (T+D ) 2 H2  - (T2+2TD+2HHq)  = 0 

A1.9 

This  is  a quadratic 

equation  in  x,  which  has  real  roots 

only  if 

its  discriminant  is  positive.  The  discriminant  is 

A = 16  d2  (T2+2TD+2HH  )2 

o 

- If  \ ( T+D ) 2 - d2l[ij  (T+D)2  H2  - (T2+2TD+2HH  )21 


A1.10 


or 

A = 16  (T+D)2  [^T2D2+^(T3+2THFq-  2TH2)  D 

+ T**  + llT2  HH  - *JT2  H2  I Al.ll 

o J 

It  Is  evident  that  A is  always  positive  for  D sufficiently  large. 
To  a first-order  approximation,  for  D large 


A = 61*  (T+D)2  T2 


Then  show 


/A  = 8 (T+D)  TD 


- / T2  + 2HH  - 2H2  \ 

V-  * ® — ) 

(T2  + 2HH  - 2H2\ 

1 + ^D ) 


A1.12 


= 8 TD2  + 12  T2D  + 8 HH  D - 8 H2D 

o 


A1.13 


With  this  approximation  and  the  approximation  D * d for  large  D, 
Equation  5 may  be  solved  to  yield,  for  large  D 


x ~ 


- 4 (T2+2TD+2HH  ) D + (8TD2+12T2D+8HH  D-8H2D) 
o — o 

16  TD 


A1.14 


This  gives 


T2  - H2 
2T 


A1.15 


and 


~ ~ D. 


A1.16 


The  second  solution  x^  Is  always  negative.  It  corresponds  to  the 
end  of  the  ellipsoid  nearer  the  radar  cutting  the  target  altitude 
plane  and  is  of  no  interest.  The  solution  x^  is  seen  to  be 
positive  for  T > H and  negative  for  T < H.  This  is  the  important 
result : 
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When  T > H,  for  any  D sufficiently  large,  the  radar  position 
will  project  inside  the  TOA  ellipse  In  the  target  altit^e  plane. 
Hence  there  will  be  exactly  one  branch  of  the  target  locus  curve 
compatible  with  large  radar  distances. 

When  T < H,  for  any  D sufficiently  large,  the  radar  position 
will  project  outside  the  TOA  ellipse  in  the  target  altitude  plane. 
Hence,  if  the  target  locus  curve  exists  at  all  for  large  radar 
distances,  it  has  two  branches.  Furthermore,  it  is  seen  that  as 
D becomes  large,  the  distance  along  the  LOP  from  the  radar  posi- 
tion to  the  TOA  ellipse  approaches  a constant  value  (Equation 
A1.15).  Also,  as  D becomes  large,  the  ellipse  locally  near  the 
LOP  becomes  more  and  more  nearly  a straight  line  perpendicular  to 
the  LOP  (cf.  Figure  A-2).  Then  for  the  radar  approximately 


2 2 

x * n - 

2T 


A1.17 


from  the  elipse  along  the  LOP  and  a sufficiently  small,  the  line 
determined  by  the  DAZ  will  intersect  the  ellipse  at  a distance 


2 2 

y » - — 2f  T'  tan  a A1.18 

from  the  LOP;  y is  approximately  constant  for  large  D.  Hence  the 
"inner  branch"  of  the  target  locus  curve  is  seen  to  approach 
asymptotically  a line  parallel  to  the  LOP  a distance  y from  the 
LOP. 

The  next  section  gives  in  a concise  form  a more  formal  and 
complex  algebraic  proof  of  the  properties  of  the  locus  lines  far 
from  OWN. 


190 


FIGURE  A-2.  DETAIL  OF  DIAGRAM  FOR  CONSTRUCTING 
TARGET  LOCUS  CURVE  . 

CONCLUSIONS 

Let 

H = height  of  TARGET 

a = differential  azimuth,  assumed  to  be  positive 
T = TOA 

The  line  parallel  to  the  line  joining  OWN  to  RADAR  and  height 
H above  it  will  be  called  the  base  line. 

Then 

1)  If  T > H the  locus  curve  has  a single  asymptote.  It  is 
parallel  to  the  base  line  at  a distance  of 

T cosa  + }/  T^  - H2  sin2a 

sina 

2)  If  H sina  < T < H the  locus  curve  has  two  asymptotes 
parallel  to  the  base  line  at  distances  of 
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T cose  ± V T2  - H2  sln2a 
sina 

3)  If  T < H sina  the  locus  curve  is  bounded. 

A.  2 JUSTIFICATION 
Let 

H = height  of  OWN  and  define  9,  r,  s,  d by  the  following 

o 


TARGET 


Here  the  points  labeled  TARGET  and  OWN  actually  represent  the 
horizontal  projections  of  TARGET  and  OWN. 

Let  x = p . 

By  the  law  of  sines 

sin  ($+ct) 
d x sina 

and 


sln$ 

s ~ x sina 


Now  T 


= + H2  + ^TTh^T)2  - fd2  + 


H 


Substitution  and  multiplication  by  T sina  results  in 
sina 


^Tx  - 


1 + (H-Hq)2  x2 


'W 


? 2 2 2 
sin^i*  + H x sin  a 


- ^sin2  (<t>+a)  + H2x2  sin2a 


A2.1 


Suppose  there  is  a not  necessarily  continuous  or  unique 
function  defined  for  small  non-negative  values  of  x,  and  taking 
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values  between  0 and  tt  — ot  (the  only  physically  realizable  bear- 


ings) such  that 
sinot 


^Tx  - ^1  + (H-Hq)2  x2)  = ^sin24>(x)  + H 


2 2,2 

x sin  a 


- ysin^Kx)  + a ^ + H2  sln2a  xc 
0 < $(x)  < it  - a 


A2.2 


If  x Is  close  to  0 


so 


slna  ~ sln$(x)  - sin ($ (x ) + a)  = - 2 cos^J>(x)  + ^ sin  j 

i ( x ) + j J 


cos  2 ~ cos(4>(x)  + 

Since  <t(x)  = - a Is  not  possible 

$(x)  = 0. 

It  follows  that  $(0)  = 0 and  that  $ Is  continuous  at  0.  By 
Taylor's  Theorem  applied  to  the  sine  function  A2.2  may  be  written 
in  the  form 


A2.4 


Vsln^a  + u ( x ) = sina  + s*lm  u(x)  + r(u(x)) 


where  R is  a function  with  the  properties 


lim 


= 0 and  R ( 0)  = 0. 


v+O 


r(u(x))  _ r(u(x))  ^ u(x) 

X ~ u(x ) * X 


Since  lim  4>(x)  = 0,  lim  u(x)  = 0 from  which  it  follows  that 
x->-0  x-*0 


p(uOoi  , 0 


llm  T77T 


x-*-0 


Note  that  u- x-^  = x ^ (2  sina  cosa  + e(x))  where  lim  e(x)  = 

X x-*-0 

It  follows  that 


v;-2 


sin^g  + u(x)  - sing  _ $ ( x ) ^co-.^  + 1 


e (x ) 


') 


2 sina 

+ • ^2  sina  cosa  + e(x)^  ^ 

_ $(x) 


|^cosa  + 6(x)^  where 


lim  6(x)  = 0. 
x->-0 


Then,  taking  difference  quotients  of  both  sides  of  A2.3, 
T sina 


')  +H2  sin2a  + (^^-)  $(x)g^$(x)^ 


$(x) 

X 


■^cosa  + 6(x)^ 


A2.5 


It  is  easy  to  see  that  if  ~p-  is  very  large,  the  right  hand  side 


of  A2.5  will  be  much  larger  than  the  left  hand  side  when  x is 


L9 


A ( jr  ^ 

small.  It  follows  that  ■ ^ Is  bounded  and  A2.5  may  be  rewritten 

T sina  * /M  sln2a~— ^X-  cosa  A2.6 

Manipulation  of  A2.6  results  in 

sina  - 2T  cosa  + ^K2-T2^  sina  = 0 

from  which 

$(x)  , T cosa  + - H2  sln2a 

x ~ sina  ' 


If  T < H sina  a contradiction  results  and  result  3 follows 
immediately . 

4>(x) 

— - — must  be  positive.  A2.7  gives  two  positive  expressions 
if  and  only  if  T2  cos2a  > T2  - H2  sin2a  or  T < H.  Note  also  that 
the  limiting  distance  from  the  base  line  is  lim  — n*  ^ x ^ = lim  - -x-^ 


x+0 


x-*0 


It  suffices  to  show  that  for  small  x 

1)  if  T > H there  Is  a value  of  $>  satisfying  A2.1 

2)  if H sina  < T < H there  is  a value  of  $ satsifying  A2.1  and 
T cosa  + | f^T2  - H2  sin2a 


$/x  < 

and  one  satisfying 


sina 


$/x  > 


T cosa  - i |/t2  - H2  sin2a 


sina 


Define  F by 


F(x,4>)  = sina ^Tx  - + (H-Fq)2x2  ) 

- ^sin24>  + H2x2  sin2a  + |4in2($+a)  + H2  x2  sin2a 

F(x, 0)  = slna^Tx  - ^ 1 + (H-HQ)2x2j-  Hx  sina  + sina  |/l  + H2x2 


s (T-H)x  sina 
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Therefore,  F(x,0)  < 0 if  T < H and  F(x,0)  > 0 if  T > H. 


Thus  if  T > H there  must  be  a t with  0 < t < j and  F(x,t) 
Claim  i)  has  been  established. 

Now  suppose  H sina  < T < H. 


Let 


Y 


T cosa 
sina 


0. 


The  second  claim  will  be  established  if  it  can  be  shown  that,  for 
small  x,  F ( x , y ) > 0. 

To  a close  enough  approximation,  sin  y = y and 

2 2 

sin  (y+a)  = sin  a +2  sina  cosa-y  (by  differential  approximation 


There fore 


x 


F(x,y)  = - sina  + Tx  sina  - x ^/T2  COSg^t  + H2  sin2a 

sin  a 


+ sina  + Tx 


cos2a 

sina 


F(x, y ) = x 


/_!_  . Jr 

I sina  y 


T2  -s„-?  + H2  sin2 


2 

sin  a 


;) 


F(x,y)  will  be  positive  if  and  only  if 
T2  j,  rp2  cos_a  + h2  sln2a 


2 

sin  a 


2 

sin  a 


if  and  only  if 
T > H sina 
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